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Abstract: We analyze the rate of convergence of the local statistics of Dyson Brownian motion to 
the GOE/GUE for short times t = o(l) with deterministic initial data V. Our main result states 
that if the density of states of V is bounded both above and away from 0 down to scales ^ t 
in a small interval of size G » t around an energy Eq, then the local statistics coincide with the 
GOE/GUE near the energy Eq after time t. Our methods are partly based on the idea of coupling 
two Dyson Brownian motions from [6], the parabolic regularity result of [15], and the eigenvalue 
rigidity results of [21]. 


1 Introduction 


Wigner ensembles consist oi N x N real symmetric or complex Hermitian matrices W whose entries 
are random variables that are independent up to the symmetry constraint W = W*. Wigner’s global 
semicircle law [28] states that in the appropriate scaling the empirical distribution of the eigenvalnes 
(Aj) converges to 

1 ^ 1 

- 2 SxAE) - PsciE) := —l{|ij|<2}V4-^^ {E e M) (1.1) 

i=l 

in the weak sense as Y ^ oo. The distribution PsdE) is referred to as the semicircle law. Wigner 
obtained his global semicircle law by computing the moments E[tr(lF"')] for each n. 

We denote by pvv(Ai,..., Aat) the joint probability density of the nnordered eigenvalues of W. If 
the entries of W are independent real or complex Ganssian random variables with variance equal to 
N~^ then the joint density is explicitly computable and is given by 

PG(Al,...,AAr) = ^]4 [Ai - Ajl^e"^^A]]IiA2/4 ^^ 2) 

^ i<j 


where /3 is 1 or 2 for the real and complex cases, respectively. Above, Zq is a normalization constant 
which can be computed explicitly. These special cases are known as the Gaussian Orthogonal and 
Gaussian Unitary ensembles (GOE and GUE). The n-point correlation fnnctions are defined by 

•••, A„) := r AAr)dAi...dAAr. (1.3) 

Eor the GOE and GUE, the n-point correlation functions have been computed explicitly by Dyson, 
Gaudin and Mehta (see, for example, [23]) nsing orthogonal polynomial techniques exploiting the 
Vandermonde determinant structure. 

Eor the case /3 = 2, the work of Dyson, Gaudin and Mehta asserts that at every fixed energy 
E e (—2, 2) in the bulk of the spectrum. 


Psc{E)' 


(n) 

iPg 


E- 


Oil 


NPsciE) ’' 


..,E 


(Xr, 


NpsciE) 


detiK{ai - aj))T=i 


(1.4) 
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where K is the sine kernel 


(1.5) 


7r(x - y) 

The rescaling by a factor of N~^ corresponds to the typical distance between consecutive eigenvalues, 
and we refer to laws under such a scaling as local statistics. There are similar but more complicated 
formulas for the GOE. 

The Wigner-Dyson-Gaudin-Mehta conjecture, or the ‘bulk universality’ conjecture, states that 
the local eigenvalue statistics of Wigner matrices are universal in the sense that they depend only 
on the symmetry class of the random matrix ensemble (i.e., real symmetric or complex Hermitian) 
but are otherwise independent of the underlying law of the matrix entries. This conjecture has been 
established for all symmetry classes in the works [7,8, 10,12,15,19]. Parallel results were obtained 
independently in various cases in [26,27]. 

We remark that when we refer to ‘bulk universality’ in this paper, we refer to the vague convergence 
of the correlation functions in the averaged energy sense of [7,8,10,12,15,19]. Results at fixed 
energy for the real symmetric case were recently obtained in [6] but we will not address this type of 
convergence in this work. 

In order to place the current work in context, we recall the three-step strategy of the proof of bulk 
universality of Wigner matrices. 

(1) Establish a local semicircle law controlling the density of eigenvalues down to the optimal scale. 

(2) Prove universality of Wigner matrices with a small Gaussian component by analyzing the con¬ 
vergence of Dyson Brownian motion to local equilibrium. 

(3) Prove universality of a general Wigner matrix by comparing its local statistics to an approxi¬ 
mating Wigner ensemble with a small Gaussian component. 

Por an overview of this strategy and a survey of recent results we refer the reader to [16]. In this 
paper we are mainly interested in Step (2). The local ergodicity of Dyson Brownian motion (DBM) 
is the intrinsic mechanism behind the universality of local statistics. 

In the present work we analyze the speed of convergence of Dyson Brownian motion for classical 
values of /3 with deterministic initial data V. Our main result is that if the density of states of V is 
bounded above and below at all scales down to .£ « t in a window of size G » -s/i around an energy Eq 
then the local statistics of Dyson Brownian motion exhibit bulk universality at time t near Eq. We 
allow for scales as small as £ = N~^, and G can shrink as long as it satisfies G » ■\/i. 

A completely analytic approach to analyzing the time to convergence of Dyson Brownian motion 
was initiated in [12], and was further developed in [13,17,19]. In these works the optimal rate t > 
was obtained when the initial data is a Wigner matrix. The key idea was to estimate the entropy 
flow of Dyson Brownian motion with respect to a global instantaneuous equilibrium state constructed 
from the invariant semicircle law of the GOE/GUE. 

Eor deterministic initial data, the study of the convergence of the local statistics of Dyson Brownian 
motion was initiated in [22]. Under some weak global conditions on the initial data, it was shown that 
the local statistics coincide with the GOE/GUE for times of order 1. In this case, the global statistics 
of Dyson Brownian motion are not close to the semicircle law and are in fact time dependent. Instead 
of comparing DBM to a global equilibrium state constructed from the GOE/GUE, an equilibrium 
state was constructed from a time dependent reference /3-ensemble specifically chosen to match the 
global eigenvalue density of the DBM. The analysis of entropy flow with respect to time dependent 
local equilibrium states was initiated in the work [29]. This method allowed for the comparison of 
the local statistics of the DBM to that of a /3-ensemble. Bulk universality for such /3-ensembles was 
achieved in the series of works [3-5] and therefore the local statistics for DBM was also obtained. 
Recently, alternative approaches to the local statistics of /3-ensembles have been presented in [25] 
and [1]. However, the key input in the proof of [22] was the universality of the local equilibrium 
measures which was only proved in [3-5]. 

As previously stated, we study the rate of convergence of DBM for short times with deterministic 
inital data satisfying only a local regularity condition. It is therefore not possible to compare the DBM 
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to a global equilibrium state as in the works [12,13,17,19,22]. One may attempt to circumvent this 
and assume regularity globally. For the class of initial data we are interested in, the local density of 
the reference /3-ensemble could have order 1 fluctuations on quite small scales and would therefore be 
quite rough. In order to complete the approach of [22], the results of [3-5] would have to be extended 
to /3-ensembles with order 1 fluctuations over small scales. While this approach may be possible, it is 
more appealing to approximate DBM locally. 

We approximate DBM locally by constructing a rescaled and shifted GOE matrix H whose eigen¬ 
values match the DBM in a small window near Eq. We then apply two ideas from previous works. 
The first is an idea of [6]: we couple the evolution of the eigenvalues of the two ensembles under the 
Dyson Brownian motion, so that their differences satisfy a system of difference equations which may 
be interpreted as a random walk in a random environment. The second idea is the parabolic regularity 
result for such systems of difference equations of [15]. In order to use the approach outlined here, two 
main ingredients are needed. These are a rigidity result for DBM with deterministic initial data after 
a short time (i.e., the analogue of Step (1) above) and a level repulsion estimate. One of our key 
observations is that the regularity of the initial data down to scales I guarantees that both rigidity and 
level repulsion occur after time t » £. 

For times of order 1, rigidity was established in the paper [21] and later refined in [22]. The 
adaptation to short times is a minor modification of the proofs there, and we will only state which 
changes are necessary in lieu of a complete proof. 

Previously, level repulsion estimates for Wigner ensembles whose entries have a smooth distribu¬ 
tion were obtained in [11]. Weaker level repulsion estimates but with no smoothness condition were 
obtained in [24,27]. The proof of [11] was modified in [6] to include the case when the Wigner ensemble 
is not smooth but instead is the sum of a (possibly non-smooth) Wigner matrix and an independent 
Gaussian part. However, these estimates degenerate as the Gaussian component becomes small, and 
would therefore only be useful for our purposes if we were interested in times of the order t = N~^. 
We establish new level repulsion estimates which show that as long as t » £ then one already has 
(essentially) as much level repulsion as one would have for times of order 1. 

Putting these ingredients together, we will prove that at a short time after the coupling is initiated, 
the eigenvalue gaps of the two ensembles coincide down to a scale with high probability. This 

proves the fixed label gap universality for DBM and also the bulk universality of the n-point correlation 
functions in the aforementioned locally averaged sense. 

As an application of our work we prove bulk universality for deformed Wigner ensembles with 
a small Wigner component. The case of a large Wigner component was proved in [22]. Another 
application of our result is the bulk universality of sparse Erdds-Renyi graphs. In a paper with J. 
Huang [20] we prove that Erdds-Renyi graphs where the probability p of each edge occuring is as small 
as p > N^/N exhibit bulk universality. The previous result obtained in [7, 9] allowed for p only as 
small as p > 

We outline the rest of the paper. In Section 2 we define our model and state our main results on 
bulk and gap universality. In Section 3 we state the local law for deformed Wigner ensembles and 
state the rigidity estimates for the eigenvalues. We also state our level repulsion estimates. Section 
4 contains the main novelty of this paper, our analysis of the Dyson Brownian motion with initial 
data V. In Section 5 prove our level repulsion estimates. In Section 6 we give the proofs of our main 
results. In Section 7 we state and prove some deterministic facts required for the proof of the local 
law, and then give a proof of the local law. We also derive the rigidity estimates from the local law. 

After completing this manuscript we learned that similar results were obtained independently 
in [14]. 

Acknowledgements. The authors thank Roland Bauerschmidt for helpful comments on a prelimi¬ 
nary draft of this manuscript. The authors thank J. Huang for pointing out an improvement in the 
assumptions of Dehnition 2.1. 
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2 Main results 


Before defining our model and stating our results, we remark that we will only state and prove our 
results in the real symmetric, i.e., /3 = 1 case. The adaptation to the complex Hermitian case, i.e., 
/3 = 2, requires only notational changes. 


2.1 Definition of model 

In this section we introduce the model under consideration. 

Definition 2.1. Let 1 = 1^ and G = Gn be two N-dependent parameters satisfying^ 

^ f <G < (2.1) 

for some ei > 0. A deterministic diagonal matrix V = diag(Vi,..., Vn) is called (£, G)-regular at Eq if 
there are constants cy > 0, Gy > 0 so that the following holds. On the interval 


Ieo,g-={Eo-G,Eo + G) 


we have 


Cy < Im my{E + ir/) < Cy 
uniformly for E e Ieo,g and i < g <10, where 

1 ^ 1 
2=1 


( 2 . 2 ) 

(2.3) 


(2.4) 


is the Stieltjes transform of V. Moreover, we assume that there is a fixed number By > 0 so that 


\VA < 


(2.5) 


for every i. 

We will use the notation 

Ie,g-={E-G,E + G). (2.6) 

Definition 2.2. The Gaussian Orthogonal Ensemble (GOE) consists of symmetric matrices W whose 
entries are independent centered Gaussians (up to the constraint wtj = Wji) with variance 

EK-] - (2.7) 

Definition 2.3. The deformed GOE consists of symmetric matrices 

Ht:=V + VfW (2.8) 

where V is a deterministic diagonal matrix and W is a GOE matrix and T > 0 is a real parameter. 

Remark. Up to a trivial constant rescaling which goes to 1 in the N ^ co limit, the solution of DBM 
with initial data V is equal to Ht at time t in law with T ~ t, for t = o(l). 

1. In a previous draft our assumption was » £ instead of G » £. The authors thank J. Huang for pointing out 
this improvement. 
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2.2 Semicircle and deformed semicircle laws 


In order to state our results on the local eigenvalue statistics we must introduce the macroscopic 
eigenvalue densities of the GOE and deformed GOE. The macroscopic eigenvalue density of the GOE 
is given by the semicircle law: 


V4 - E2 

Psc{E) := 1{\E\<2} - ^ -■ 


(2.9) 


While the macroscopic eigenvalue density of Wigner (and generalized Wigner) matrices also follows 
the semicircle law, the deformed GOE follows a deformation of p^c, the so-called free convolution of 
the semicircle law and the empirical measure of V. We define it through its Stieltjes transform. We 
let be the solution to 


m 


(N) 

fc,T 


N 




Vi- z- Tm\^l,{z) 


fc,rv 


( 2 . 10 ) 


Above, Vi denote the entries of V. The properties of the above equation are well-studied. It is known 
that there is a unique solution to the above equation and that is the Stieltjes transform of a 

measure which has a density p[^^- This density is compactly supported and analytic on the interior 
of its support. We refer the reader to, for example, [2] for details. We emphasize that and p^^^ 

and their qualitative properties depend on N. For example, the density can become quite rough 
as —> 00 . For notational convenience we suppress the superscript and denote 


mic,T{z) := Pic,T{E) := 


( 2 . 11 ) 


When we wish to emphasize the ^dependence we will use rn^^^ and instead. 

To state our result on gap universality we define the classical eigenvalue locations of the GOE and 
deformed GOE by 

i r'yT,i i rui 

^ = J_^ ftc,T(i^)dE, - = psc(^)dE. (2.12) 

The classical eigenvalue locations of the deformed law depend on N but we again suppress this in our 
notation. 


2.3 Bulk and gap universality 

For a deformed GOE matrix H = V + 'JTW we denote the re-point correlation function by . It is 
dehned by 

p5?^(Ai,...,A„) = f pi^^)(Ai,....,A,v)dA„+i...dA,v (2.13) 

where is the joint density of the unordered eigenvalues of Ht- The re-point correlation functions 

(n) 

of the GOE are denoted by Pqq^ and are defined similarly. For the joint eigenvalue density of the 
GOE we have the explicit expression 

= (V) (2-14) 

^GOE i<j 

where is a normalization constant. Our main result on bulk universality is the following. 

Theorem 2.4. Let 

Ht = V + VtW (2.15) 

be a deformed GOE matrix. Suppose that V is {I, G)-regular at E and that N~^G > T > for 
some e > 0. Let O e C'“(M"') be a test funetion. Fix a parameter b = N^/N for any c > 0 satisfying 
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c < e/2. We have, 


nE-\-h n 

lim 0{ai,...,an) 




(n) 

-Pt 


E' 


CX\ 


Np\^_l{E) 

7-1// , 


..,E' 


CXr 


ipsc{E")y 


;Pgoe {+ 


Npsc{E") 


,...,E" + 


NpsciE") 


dE' 

dai...dan-^ = 0 (2.16) 


for any E" e (—2, 2). 


Remark. The upper bound on the size of the averaging window can be removed as long as it is 
contained in Tg c /2 a-iid one replaces pic^j-^E) by pic,TiE'). This is due to the fact that the macroscopic 
density /9fc,r varies on the scale T~^; see Lemma 3.2 below. 

Remark. The scaling factor Pf^^{E) satisfies c < p^^^{E) < C for all large N] see Lemma 3.2. 

For the gap universality we have the following. 


Theorem 2.5. Let 

Ht = V + VtW (2.17) 

be a deformed GOE matrix. Suppose that V is {£, G)-regular at E and that GN~^ > T > N^£ for some 
e > 0. There is a > 0 such that the following holds. Let O e (//“(M”) he a test function. Let i be 
an index so that the ith classical eigenvalue of Ht (defined in (2.12)) lies inX^Qj^. Let ii,...,in e N 
with ik < for each k. Then for N large enough we have 




O - A,+iJ,...,iVp;5(7T,)(Ai - A,+ij)^ 

- [O {Npsc{pi){\i - Xi+ifis,Npsc{pi){\i - Aj+i„))] 




(2.18) 


2.4 Bulk universality for deformed Wigner ensembles 

In this section we summarize how Theorem 2.4 can be used to prove bulk universality for deformed 
Wigner ensembles with a small Wigner component. Deformed Wigner ensembles are of the form 

Ut-.= V + VtW' (2.19) 

where F is a deterministic diagonal matrix and W' is a Wigner matrix. In [22], bulk universality was 
proven for such ensembles for T = 1 by following the three step strategy outlined in Section 1. Using 
our result one can use the three step strategy to show that if the density of states of V is bounded 
above and away from 0 down to scales £« T in a window of size G » \/T around Eq, then Ut 
exhibits bulk universality in the averaged n-point correlation sense at the energy Eq. As previously 
stated, we have shown how to adapt the proof of Step (1) in [21,22] to short times in this paper. The 
adaptation of Step (3) to deformed ensembles was also achieved in [22], and the adaptation to small 
times is straightforward. We omit it in this paper. This completes the proof of bulk universality for 
such ensembles. 

The above proof can clearly be extended to the case when V in (2.19) is random and independent of 
W, as long as the regularity properties are uniformly satisfied with high probability and the empirical 
Stieltjes transforms of V are close to some deterministic Stieltjes transform. This is a further simple 
modification of the arguments of [22] and we omit the details. 


3 Local deformed semicircle law, rigidity and level repulsion 

3.1 Conventions and notations 

In this subsection we introduce some notation and conventions which will be used throughout the 
paper. We introduce two parameters f = f,N and p = We let 

ao < ^ < AologlogV, p = (log(iV))‘^i (3.1) 

for some constants oq > 2, Aq > 10 and Gi > 1. 
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Definition 3.1. We say that an event has {^,v)-high probability, if 


o>0 (3.2) 

for all large N. We say that an event Hi holds with {f,,v)-high probability on an event H 2 if 

P[H^ n H 2 ] < (3.3) 


3.2 Deformed semicircle law 

For reasons that will be clear in Section 4 we modify the definition of Ht slightly and introduce 

Ht ■= + {1 - e-^)^/^W. (3.4) 


The solution of DBM with initial data V is equal to Ht in law for each t. Up to a trivial rescaling 
Ht corresponds to the deformed model Ht with T ~ t, for times t = o(l). Note that the assumptions 
(2.3) on V are essentially invariant under rescaling V for t = o(l). 

The deformed semicircle law satisfies the equation 


m 


(TV) 

fc,t 


(z) 


1 ^ 1 

ly_ t _. 

^ iWi e-^/^V) - z - {1- {z) 


(3.5) 


As above, is the Stieltjes transform of a measure that has a density that is compactly 

supported and analytic on the interior of its support. For notational convenience we will suppress the 
superscript and use the notation 


mic,t{z) ■= Pic,t{E) := (E). 


(3.6) 


When we wish to emphasize the A'-dependence we will use and p^^^ instead. 

We now collect some properties of the density pfc^t and We will assume that V is (£, G)-regular 

at Eq. For w > 0 satisfying 

u; < ei/10 (3.7) 

where £i is the constant appearing in (2.1) we define the time domain 


% = {t:m‘^ <t< N-‘^G}. 


(3.8) 


We defer the proof of the following lemma to Lemma 7.2 in Section 7. 

Lemma 3.2. Fix w > 0 and 1 > q > 0. There are constants c,C depending only the constants in 
(2.3) such that 

c < ptc,t{E) < C, \pic,tiE)\ < j (3.9) 

for all E e XEQ,qG-,t ^ N large enough (where large enough depends on uj, q and the constants 

in (2.3)). 


3.3 Local law and rigidity 

For any L > 12^ and 0 < g < 1 we define the spectral domain 


VL,q := Vi u V 2 (3.10) 

where 

Vi := {z = E + ip : E e lEo,qG, <Np< lOA} 

V 2 := {z = E+ ip : |F;| < lOA < Np < ivio^v+2j^ 
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The domain T>i is where we prove the local deformed law down to the optimal scale rj > N~^, using 
our assumptions on V. On the domain P 2 we are only concerned with the behaviour of ttiat when z 
is either far from the real axis or far from the spectrum of Ht- 
We define the empirical Stieltjes transform of Ht as 


'mN,t{z) = —tr 



(3.12) 


Theorem 3.3. Let Ht he as in (3.4) where V is i-regular on I. Let 




^0 + 0(1) 
2 


log log N. 


(3.13) 


Fix 0 < q < 1. There are constants ly > 0 and ci depending on Xy, the constants appearing in (2.3), 
Aq and Cl in (3.2) and the choice of q, so that the following holds with L > 40^. For any z e 
and t e Tco 

\mN,t{z) - m[^j{z)\ < (3.14) 

with {f^,u)-high probability for N large enough (here, large enough depends on the choice of v, oj and 

q)- 

The proof follows closely that of the same result in [21,22]. We summarize the changes in Section 

7. 

From this local law we derive the following result on the rigidity of the eigenvalues of Ht. We 
define the classical eigenvalue location jt^t by 


i 

N 


rji.t 

J—00 


Pfc,t(u)du. 


For a fixed constant q > 0 satisfying 

0 < g < 1 


and a time t we define the bulk index set by 


(3.15) 

(3.16) 


•^q,t ■ {i • 'Jitt ^ (3T7) 

We abbreviate 

;= K{Ht) (3.18) 

where Aj(T) denotes the ith largest eigenvalue of a matrix A. The following lemma shows that the 
definition of Aq^t is essentially independent of the chosen time t. 

Lemma 3.4. Let 0 < < ^2 < 1 o-iT'd a; > 0. Then for N large enough, we have for all s,t e 7Ii, 

Aq,,s C Aq^^f (3.19) 

We defer the proof to Section 7. 

Theorem 3.5. Fix 1 > q > 0. There are constants y > 0 and C 4 > 0 depending on Xy, AQ,aQ, Ci 
and the choice of q such that the following holds. For i e Aq^t we have 




(3.20) 


with (^, y)-high probability for all t e Tu,- 

Deducing Theorem 3.5 from Theorem 3.3 is a minor modification of the usual approach given in, 
for example, [8], and so we will not provide all of the details. This is sketched in Section 7.3. 




3.4 Level repulsion estimates 

The following level repulsion estimates will be an important tool for our proof. The proof is given in 
Section 5. 

Theorem 3.6. Let <5 > 0 and r > 0 and let denote the eigenvalues of the deformed GOE matrix 
Ht with V (£, G)-regular at Eq. Eix cv > 0 and 1 > (; > 0. Eor t and i e Ag^t we have: 

P[|Ai,t - Ai+i,t| < e/N] < (3.21) 


4 Optimal speed of convergence of Dyson Brownian motion 


4.1 The approximating GOE ensemble and statement of result 


In this section we implement our main strategy in proving gap universality with a fixed label by 
approximating Ht by a GOE matrix around a fixed eigenvalue. Let IE be a GOE matrix. We denote 
the semicircle law by 

p^c{x) := (4.1) 

and the semicircle law for alE + 6 by 


Psc \x) : l{|a;- 6 |< 2 a} 2a?'K 


— {x — by 


The classical eigenvalues are defined by 


(4.2) 


A 

N 


l: 




— 2 ( 2+6 


rsc 


{u)du 


(4.3) 


and we abbreviate pk = ■ The following lemma shows that at a fixed time we can locally 

approximate the classical eigenvalues of Ht by those of a GOE matrix. We omit the proof as it is 
elementary. The only ingredients are the bounds of Lemma 3.2. 


Lemma 4.1. Let ko e Ag^t be a bulk index, and t e Eix a > 0 and let jo satisfy aN < jo < 
(1 — a)N. Let a and b be 

6 = (4.4) 


a = 


Pfc,ti'yko,t) 

Then there is a constant C depending only on a and the constants in (2.3) so that 


C 


\Pjfl - 7i+fco,d ^ for\j\ < \fNi 


(4.5) 


In the following we will implement a coupling argument to match the eigenvalue gaps of Ht by 
those of a GOE matrix down to a scale We define A,- as the solutions to the SDEs 


dXJt) = \ —dBi + I — V --Ai(t) 1 dt (4.6) 

with initial data 

Ai(0) = Ai(E) (4.7) 

the ordered eigenvalues of V. Eor each fixed time t the vector A is distributed as X{Ht). Fixing a time 
to > 0 we define fii{t) for t > to as the solution to the SDE 


dfLi{t) 



fly_i_ 

-Pk{t) 


1 


Pi{f) dt 


(4.8) 
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with initial data /ii(to) = ^iiaW + b) with W a GOE matrix and constants a and b which will be 
chosen according to Lemma 4.1. The Brownian motions appearing above are the same as in equation 
(4.6). This idea of coupling two DBMs originated in [6] and is crucial to our work. 

For fixed times t the vector p, is distributed as the eigenvalues of atW + bt with deterministic 
constants at and bt which are easily computed and satisfy c < at < C and |6t| < C. 

The goal of the remainder of the section is the proof of the following theorem on the differences of 
the gaps (Aj - Xj){t) - - i^j+jo-ko){t) for a time t = to + N^-^. 


Theorem 4.2. There exists universal constants a > 0 and ci > 0 such that the following holds. Let 
V be {i,G)-regular at Eq. Let A(t) be the solution to the SDE (4-6) with initial data V as described 
above. Fix uj > 0, 1 > q > 0 and let to = N'^°/N e %j. Fix a bulk index ko e Agg^, a > 0 and choose 
a,b,jo in Lemma 4-1 for t = to, and let W be a GOE matrix independent of Htg. Let p be the 
solution of (4-8) with initial data aW + b at time to as described above. 

Let 0 < uj' < To. Then, there is an event E with probability > 1 — on which the 

following holds. 


{Xj+koil'O + I') ^j'+fco(to + t) {pjj,.jg(to + t) + t))| ^ pjauj'+l ’ 

_ ]<guj'/50 jyaj'/lO 


for \j\ + \f\ < iV‘"'/4000^ ^ 


N 


N 


(4.9) 


Over the next few subsections we will establish the estimates we will need to prove Theorem 4.2. 
In the last subsection we will provide the proof. The strategy is roughly as follows. Due to the fact 
that we have coupled the DBM of both ensembles, their differences satisfy a system of difference 
equations. The Holder regularity result for such parabolic systems of [15] states that such solutions 
(neglecting edge effects) v satisfy an estimate of the form \vi{t) —Vi+i{t)\ < A^“^||u(0)||oo after a short 
time t. If we knew that the classical eigenvalue locations of our two ensembles matched throughout the 
spectrum then by the optimal rigidity results we would have ||u(0)||oo < with overwhelming 

probability, and we would deduce that the gaps of the two ensembles coincide. However, as we only 
have the matching for roughly eigenvalues we will have to cut off the non-matching eigenvalues. 
While this cut-off disrupts the location of the eigenvalues, the gaps remain unchanged down to a scale 
However, due to the cut-off, the differences of the two eigenvalues will only be an approximate 
solution to the system of difference equations that we have the parabolic regularity result for. The 
error term will be controlled by the finite speed of propagation estimates obtained in [15], and we 
will see that the parabolic regularity result in fact applies to our approximate solution as well. We 
therefore obtain the desired statement about the gaps. 


4.2 Rescaling, time shift and regularization 

It is natural to rescale the problem and work in the microscopic scaling. We remark that we only work 
in microscopic coordinates for the proof of Theorem 4.2. In all other sections of this paper we state 
and prove our results in the macroscopic scaling in which the typical distance between consecutive 
eigenvalues is 0(1/A^). For this proof we also introduce a time shift and a relabeling which sets 
to = ko = jo = 0. To be more precise, we define the variables 

Xj{t) = NXk^+j{to + N~^t), yj{t) = Npj^+j{to + N~^t). (4.10) 

They satisfy the stochastic differential equations 

- V2dB, + 2 (4.11) 

d»(t) . V2dB, + 2 (4.12) 

where are independent Brownian motions. Note that the for the x^’s and y^’s are the same. 

We first show the existence of a set Qs of good paths on which rigidity holds for the eigenvalues 
Xk and Vk. 
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Lemma 4.3. Let u' he as in the statement of Theorem 4- For any <5 > 0 there is an event Qs with 
P[^ 5 ] > 1 — for all large N such that the following holds. 

c\i-j\ < \xi{t) - Xj{t)\ < C\i-j\, \i-j\ > N‘^'\\i\ + \j\ < t < 

(4.13) 

\xi{t) - Xj{t)\ < GN^'^\i-j\, \i\ + \j\ < GN^-‘^/^°,0< t < 

(4.14) 

and also the same estimates with the Xj ’s replaced by yj’s. The constants c, G depend only on the 
constants appearing in (2.3). Also on Qs 

Im {mN{E + ir])] < G, E e Tg^c, 5 > r] > N^^'^/N (4.15) 

and the same for the Stieltjes transform of the matching GOE ensemble. 

Proof. The result for the xj's for fixed times t on a set of high probability follows from the rigidity 
estimates of Theorem 3.5. The extension to all times is a minor stochastic continuity argument almost 
identical to the proof of Lemma 9.3 in [15], with the input of the level repulsion bounds of Theorem 
3.6. The result for the y^-’s is similar. □ 

We will also need to regularize the dynamics.We define the regularized eigenvalues as the solution 
to the equations 

dxk = V2dBk + 2 — , . _ -dt - ■:^Xk{t)dt, Xk{0) = Xk{0), (4.16) 

XjV') ~r Skj zjV 

The constant Skj is defined as ekj = £ for k > j and £kj = —£ for k < j. Here s = We make 

the same definition for the y^’s and introduce the regularized y^’s. 

The following lemma is based on [6]. The main inputs are the level repulsion estimates stated in 
Theorem 3.6. 


Lemma 4.4. Choose e as 


e = N 


-500 


Then there is an event Ei with probability > 1 — on which 

\xk{t) - Xk{t)\ < 

for every k with k + ko e Ag^tg ond all 0 < t < G‘^N^~^I‘^. 

Proof. The rescaled difference qi{f) := (^Xi{t) — Xi{t)) satisfies the equation 


(4.17) 

(4.18) 


dqi = Qi{t)dt, with Lli{t) = ^ 




jjti 


{Xi{t) - Xj{t)){Xi{t) - Xj{t) + Sij) ' 


We bound 


IE sup 

0<t<G2Ari-W2 


r Lli{s)di 

Jo 


< 


sup 


sup ( E 


0<t<G2Ari““/2 


1 \ 


2/3 


I CC i X j 


|3/2j 


E. 


\Xj — Xn + e. 




By (3.21), as we have restricted our attention to i + k^ e Ag^t, we have 

1 \ 2/3 


E- 


\Xi Xi.^-1 


13/2 


< N 


(4.19) 


1/3 


(4.20) 


(4.21) 


and for any r > 0 
E 


\xi+i -Xi + eP 


1/3 


< E 


-2-r 


\xi+i - Xi + e|2 


\ 1/3 

j <iVe(2-0/3|iog^|i/3. 


The claim then follows from the Markov inequality and our choice of e. 

For the y^’s we have even better level repulsion and so Lemma 4.4 holds for the y^’s as well. 


(4.22) 

□ 
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4.3 Cut-off of non-matching eigenvalues 

Let now w' > 0 be as in the statement of Theorem 4.2 and let 

K = 

Define Xj for \j\ < K by 


dxj{t) = V2dBj + ^ 


1 


1 


lTtk,\l\<K 


Xj{t) - Xi{t) + Eji 2N 


dt - —Xj{t)dt, Xj{0) = Xj{0), 


(4.23) 


(4.24) 


and a similar definition for yj. We prove the following lemma to control the effect of the cut-off. The 
main inputs are the level repulsion estimates of Theorem 3.6 and the rigidity estimates of Lemma 4.3. 

Lemma 4.5. Let 1/2 > 5 > 0, and ti < There is an event J- 2 ,& with P[J- 2 ,( 5 ] > 1 — 

on which the following estimates hold. We have, 


sup \Xa{t) - Xh{t) - {Xait) - Xb(t))| < 


K\ti + l)\a-h\ 


0<t<ti 


{K-\a\ + l)V2(i^ _ |5| + 1 ) 1/2 


(4.25) 


for all |a| < |6| < K, and then same the estimates changing Xj for xj, and also the same estimates 
for the yj’s. 


Proof. By Lemmas 4.3 and 4.4 we can assume that Qs and J-i hold. Define for |f| < K, 

Wi{t) = e^l^‘^^\xa{t) - Xb{t) - {Xa{t) - Xfe(t))) 

to be the rescaled difference. We have, 


(4.26) 


(t) = fe^/(2^) 2 


aia('S) - Xbis) 


{Xa{s) - Xj{s) -h eaj){Xb{s) - Xj{s) + Ebj) 


ds 


E 


Xa{s) - Xbis) 


K+K^>\j\>K 


(Xais) - Xj{s) + Eaj)iXbis) - Xj{s) + Ebj) 


ds 


f 


,s/i2N) 


E 


Xais) - Xbis) 


\j\>K+K^ 


iXais) - Xjis) + Eaj)iXbis) - Xj(s) -h Ebj) 


ds 


= ^1+^2 


(4.27) 


We have split the sum into the terms involving j near the edge K and those with j away from the 
edge K. We will use rigidity to estimate the terms away from K and the level repulsion estimates for 
those near K. For A 2 we get, using Cauchy-Schwartz, 


A 2 I < CK^\a- 


oij'f E 

Vli|>K+K« 


1/2 


1/2 


1 


\Xjis) - Xa(s)|2 


E 


1 


W\>k+k^ 


|xj(s) - a:fe(s)|2 


ds. (4.28) 


We have the estimate 

1 


E 


< 


JL-|a| + l 


/v2(5 

Im [mNixais)+iK^^'^iK-\a\ + l)/N)] < —— . (4.29) 


JL — |a| -I- 1 


Hence, 

We estimate Ai by 


l^2| < 


K^^\a - b\iti + 1) 


iK - |a| -h l)^/‘^iK -\b\ -h 1)1/2' 


l^il ^ 


If 

1 Jo 


iL-|a|-hlJo (xk(s) - Xi^+i(s) - e)2 ix-xis) - X-K-iis) + e)"^ 


ds 


(4.30) 


(4.31) 
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(4.32) 


We have estimated the terms ||j| — K\ < using, for j near K, 

la^aCs) - Xj{s) + Eajl > K~^{K - |a| + l)|a;x(s) - a;x+i(s) - e 


and also 

|X6(S) - Xj{s) + £bj\ > \xk{s) - 
and then similar estimates for j near —K. We have 


[Jo {xk{s) - XK+i{s) - e)2 

by the level repulsion estimate (3.21) and the same for —K. By the Markov inequality there is an 
event of probability at least 1 — K~^ on which it less than . This yields the estimate for the 

Xfc’s. Since the event J-i of Lemma 4.4 holds we also get the estimate for differences of the x^s and 
Xfc’s. The same proof works for the y^’s. □ 

Let Ua{t) = e^/^‘^^\xa{t) — ya{t)), |a| < K. For u we have the equation 

duj{t) = ^ Bji{t){xj{t) - xi{t) - {yj{t) - yi{t)))dt 

= Yi + Cjit)dt (4.35) 


K+i{s) - e| 


< tiK^ 


(4.33) 


(4.34) 


where 


Cj{t) = Y Bji{t){{xj{t) - xi{t)) - {xj{t) - xi{t)) - {yj{t) - yi{t)) - {yj{t) - yi{t)) (4.36) 

l^j,\l\<K 


and 


{Xj{t) - Xi{t) + £ji){yj{t) - yi{t) + £ji) ■ 

(4.37) 

Let Vj, \j\ < K he the solution to 


^Vj{t) = ^ Bji[f){vj{t)-vi{t))dt 

(4.38) 



with initial data 


o' 

II 

0^ 

(4.39) 

By the Duhamel formula we have 


Ua{t)-Va{t)= r Y ^ap\s,t)^p{s)ds 
\P\<K 

/ 7 -.\ . . 

(4.40) 


where is the semigroup associated to the equation (4.38), and denote its matrix elements. 
To estimate the RHS we require the following two lemmas. 


Lemma 4.6. Let (5 > 0 and 0 < ti < GN^ There is an event J- 3^5 with probability at least 

P[J- 3 ^ 5 ] > 1 — so that the following estimates hold. 


lea(i)| < 


iL — |a| + 1 


(|-6a,a+l(i)| + |.Ba,a-l(i)|) 


(4.41) 


for every |a| < iL and 0 <t <ti. 
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Proof. Let 5 > 0. We can assume that the events ^ 5 , T\ and of Lemmas 4.3, 4.4 and 4.5 hold. 
Proceeding directly from the definition, we estimate the numerator using Lemma 4.5 and obtain 


|Ca(t)| <2(l+ti) 




\a-j\ 




< 2(1 + + \Ba^a-l\) ^ 








Lf — |a| + 1 

1 


-K 


2S 


{K - |a| + 1)1/2 |« _ j\^K -\j\ + 1)1/2 


(4.42) 


where we have used Lemma 4.3 in the last step to estimate the terms |a — j\ > K^, \ j\ < K. We now 
estimate the sum 

1 ^ 1 


E 


< 


E 


IjI+ 1)^'^^ |j|<i^^X-|a| + l>2(i^-|i| + l) 1“ lil + 1)^'^^ 




1 


Clearly, 


\j\<K,K-\a\ + l<2{K-\j\ + l),j^a 

c 


A 2 < 


a-jl(K-ljl + iy /2 

log(iL). 


=:Ai + A2. (4.43) 


{K -\a\ + 1)1/2 

For Ai note that if {K — |a| + l)/2 > {K — \j\ + 1) then \j\ — |a| > {K — |a| + l)/2 and so 

C ^ 1 


(4.44) 


4li < 


E 




< 


c 


E 


< 


ciog(ii:) 


(K - l«l + 1)'/' ^ ^ (^ - l«l + 1)'/' ■ 


This yields the claim. 


(4.45) 

□ 


4.4 Finite speed of propogation and Holder regularity 

We require the following finite speed of propogation estimate, the proof of which is a minor modification 
of the proof of Lemma 9.6 of [15]. 

Lemma 4.7. Fix 5i > 0 and ^2 > 0 and 0 < ti < . Suppose that 

k-a 


Bij > 


h -2l 


and 


sup 

0<M<K 


f ' i 5i7(s)ds < 

l + ti L M ^ 


hold. Then for any 0 < s < t < ti we have the estimate 


{SB) - \p-a\ + l 


(4.46) 

(4.47) 

(4.48) 


We can now estimate the difference between u{t) and v{t). Again, we will rely crucially on the 
level repulsion estimates of Theorem 3.6. 

Lemma 4.8. Let 5 > 0 and 0 < ti < . For every joj < there is an event F 4 ^s,a with 

P[F 4 , 5 ,„] > 1 - on which 


, /X /M (1+ 

sup \Uait) - Ua(t)| < --. 


(4.49) 
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Proof. In order to apply Lemma 4.7, we must check that the conditions (4.46)-(4.47) are satisfied 
on a set of high probability. Without the supremum, (4.47) with exponent d > 0 holds on a set of 
probability at least 1 — for fixed M. The supremum can be replaced by a dyadic choice of 

M = 2^ for k < ClogN. Hence we have that (4.47) holds with exponent 5 > 0 on a set of probability 
at least 1 — The condition (4.46) is a consequence of rigidity (Lemma 4.3) which holds on a 

set of much higher probability. Let be the event that all of this holds as well as the event of 
Lemma 4.6. Then P[J7i] > 1 — . From (4.40) we have, for all times 0 < t < H, 

\^S^Hs,t)\\^p{s)\ds. (4.50) 

Using the estimates that hold on provided by the previous two lemmas we get 



E 




2 f 

Jo 


\P\<K 


< 


2(1 + ti)5/2i^45+l/2 E , sup 


lE[JJp,p+i(0] 


\p\<K+l 


o<t<ti (|a -p\ + 1 ){K - \p\ + 1) 


< (1 + 


K5S 

1 ^' 


(4.51) 


Hence, there is an event of probability at least 1 — K~^ on which the RHS of (4.50) is less than 
This yields the claim. □ 

We now state the Holder regularity theorem (Theorem 10.3 of [15]) from which Theorem 4.2 will 
follow. In order to do so, we need to introduce the following two conditions. 

(Cl)p^ At time a we have 


sup sup -j— 

0<s<cr' 1<M<K 1 + |s — cr 


if ,S 


u < 


\i\<K\j\<K 


for every a' e {crH + a} where 

H = {-2"”*(1 + 2"^) :0<k<m< ClogN} 
(C 2 )p 2 For every 0 < s < u we have 


(4.52) 


(4.53) 


K-P 2 

Bij{s) > - -(4.54) 

and 

< Bijis) < for |f - j| > C'NP^ (4.55) 

Theorem 4.9. There is a universal constant o > 0 such that the following holds. Let v{t) be a solution 
to the equation (4-38). Let C 3 > 0 be given, and let a e \K'^^, Suppose that (C'l)p^ and {C 2 )p^ 

hold for Pi and p 2 small enough depending on C3. Then we have 

sup |uj(t) - Ujv(t)| < C'cj”“||u(0)||oo (4.56) 

|j|+|i'|<o-2/3 


for all t B [a — cr]. 

4.5 Proof of Theorem 4.2 

Recall our choice of iL = Let ti = g^j^d 6 = 1/10. By Lemmas 4.5 and 4.8 there is an 

event with probability > 1 — on which 

sup sup \xj{t) - xy{t) - {yj{t) - yy{t)) - {vj{t) - Vf{t))\ < - 7 ^, (4.57) 

|i| + |i'|<Xi/2000 0<t<ti B / 
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where Vj was defined in (4.38). We now apply Theorem 4.9 to the difference Vj{t) — Vji{t). We take 
C 3 = 1/10. By Lemma 4.3 we know that {C2)p^ occurs on a set of high probability. The other condition 
(C'l)pj occurs on a set of probability at least 1 — If there were no double supremum this 

would follow immediately from the level repulsion bounds (3.21) and Markov’s inequality. The double 
supremum can be replaced by a dyadic choice of s = 2““cj and M = 2^ for integers a,b < C log N. 
Hence we see that there is a constant ci s.t. on an event of probability at least 1 — we have 

sup |uj(t) - Uj/(t)| < CiL”"/^°||t>(0)||oo. (4.58) 

|i| + |T|<i^l/30 

for t e By the rigidity estimates for and HtQ and our choice of a and h 

(recall Lemma 4.1) we have that ||u(0)||oo < on a set of high probability. This concludes the 

proof of the theorem. □ 

5 Level repulsion estimates: proof of Theorem 3.6 

In this section we prove Theorem 3.6. Previous level repulsion estimates for Wigner ensembles were 
obtained in [11]. These estimates required a smoothness condition for the law of the entries of W . The 
dependence of the estimates on the smoothness of the matrix entries was analyzed in the work [ 6 ] . As 
a result, the smoothness condition could be relaxed at the cost of introducing an extra error prefactor 
in the level repulsion estimates depending on the size of the smoothness component of the matrix 
entries. 

Due to the smoothness condition on the matrix entries, the estimates obtained in [6,11] are not 
sufficient for our purposes. If we were interested only in times t = N~^ for small e > 0, then these 
estimates would suffice. As we are interested in much smaller times, we must obtain new estimates. 
Our overall strategy is similar to the method of [11] but we will take advantage of the fact that the 
Wigner part of our ensemble is a GOE matrix which allows us to explicitly compute the expectations 
of Green’s function elements that arise in the proof. 

We remark that a form of level repulsion was obtained in [24, 27] for Wigner ensembles under 
no smoothness assumptions on the matrix elements. However, the estimates obtained are not strong 
enough for our methods. As an additional remark, we note that gap universality implies a weak form 
of level repulsion; still, such an estimate resulting from gap universality would not be enough for our 
proof. 

The estimates we will obtain are not completely optimal but will suffice for our purposes. In this 
section, it will be more notationally convenient to work with the definition 

Ht = V + VfW. (5.1) 

Before we proceed we remark that we have switched back to the macroscoping scaling, so that the 
typical distance between eigenvalues is 0{1/N), instead of the microscoping scaling which was only 
used in the proof of Theorem 4.2. 

Theorem 3.6 will be derived as a consequence of the following theorem. For a 1/2 > e > 0 and 
E 6 lEo,qG let 

Ie:=[E- e/N, E + e/N] (5.2) 

Theorem 5.1. For k = 1,2, and 1 > g > 0 and all 6 > 0,r > 0 we have for E e 

F[Nx, >k]< ^5 3 ) 

Proof. Our starting point is the inequality 

F[Nx^ > 1] <£2E[(Im [mN{E + i 7 j)]f] < ^Y,n\Gn\\G,,\] < ^ ^ ] V2 (5.4) 

ij iJ 

where Nrj = £. We will compute the expectations appearing on the right. 
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We have the formula 


Gii = i Vi + VTwii - E -ir]-^ djij - 


(5.5) 


where 


and 


o’ ^3 




N{xf - E) 


N‘^{\y - EY + £2 - EY + £2 


^ I Xj 


U) 


(5.6) 


(5.7) 


where hik = Vi5ik + VTwik are the entries of Ht, is the jth normalized eigenvector of the ith 
(i) 

minor of EIt and X)' is the jth eigenvalue of the ith minor. 


From the formula 


i r°° ■ 

h““du, Im la] > 0 

« Jo 


we have 


rCO pi 

\Gu\^ = 

Jo Jo 


exp 


-(u + ^^) 2 Cj^j + i(u — v)^ dj^j + i{u — v){E — Vi — VTwu) — (n + v)r] 


(5.8) 


dndu. 


(5.9) 


We denote by E, the expectation over the ith row of W. This is a Gaussian integral which we can 
compute explicitly, as conditionally on h!^\ the variables xj = ^ are independent 

Gaussian random variables with variance T. Hence, 


J 'OO roo r r N—1 ^2 

exp - V -^(1 + 2 T(tt + u)cj + i(a - u) 2 rdj) 
0 Jo Jr^ L ^ 


E, 


— tY + i(u — v){E — Vi — \Y2TN ^/^y) — {u + v)ri 


f'CC pec ^-{u-v)'^ ^-{u+v)7]+i{u-v){E-Vi) 


J 'OO rCC 

0 Jo n« 


(1 + 2T{u + v)ci + i2T{u — v)diYG 


dydx 
_ (27r)^/2 

dttdu. 


dndu 


(5.10) 


Note that we also integrated out the variable wu. After the change of variables 2T{u+v, u—v) —> (n, v) 
this becomes 

rco pco ^-v'^-^^-ur]{2T)~'-+\v{E-Vi){2T)~'^ 

(5.11) 


1 

E*[|Giip|] = ^ daj du 
Let 1 > qi > q. We let Qs^i be the event 

Q5,i '■= 


W 1 Y + UC 1 + iudi)V2 


{Z:X(77«)eX,} 


< 


(5.12) 


We know that Q^^i holds with (^, i/)-high probability. Furthermore, Qs^i is independent of the ith row 
and column of Ht- Let us first get a bound for E[lc;|. On this event we know that there are 

at least 4 eigenvalues in and so there are at least 4 indices aj s.t. Cq. > s~^. Ghoosing such indices 
we bound 


Glgc 


rp2 

< \g^ e 


■ 

Jo nj=iii 


1 


du X 


+ Cq,.u|V2 J_oo 


r 




by our choice of Cq, .. Hence, 


< eN^< e. 


(5.13) 

(5.14) 
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(5.15) 


We now work on the event Our goal is to obtain 


mgsJGu\‘^] ^ 


1 




{Vi-E)^ + T^' 

This will be the result of obtaining two bounds on the integral on the RHS of (5.11) and taking the 
minimum at the end. We begin the derivation of the first bound. We choose distinct indices a* and 
/3i and use the bound 

1 


1 


1 + CaiU|V2|l + Ca,2u|V2|l + c, 


0:3 


U\ 


r 


du- 


1 


(1 + |d/3,u|)V2(l + \dp^v\YI^{l + \di3^v\Y ■ 

We now define a random variable Ai as follows. On the complement of it is 0. If holds and 
there are at least 20 eigenvalues of greater than E + e/N we set it equal to A|A;(//y^) — E\ 
where I is the index of the 20 th such eigenvalue (i.e., the smallest I s.t. Xi- 2 o{H^^) is still greater 
than E + e/N). If 0s,i holds and there are less than 20 eigenvalues greater than E + e/N then there 
are at least 20 eigenvalues less than E — e/N and we set it equal to N\\i{hIY'^) — E\ where I is the 
index of the 20th such eigenvalue. If Qsi holds and Ai < 1 we instead redefine it so Ai = 1 (this is 
just so M —>■ AY is increasing). By rigidity, the decay of the entries of W and our assumption that 
||F|| < N^^ it is easy to see that 

E[Af ] < N^ (5.17) 

for any M > 1 and 5 > 0. By dehnition of Qs,i we can always find distinct indices so that > e/(2A^) 
and > (2Ai)”^ (for example, take a* so that Aq,. are the closest eigenvalues to E and then the 
Pi’s to be the next eigenvalues that are closest to E but outside the interval Xg). With this choice we 
get from (5.16) that 

(5.18) 


mgs,\Giin ^ 

For the other bound we start with the expression on the RHS of (5.11) and integrate by parts in v 
twice using 


- 2 ri d 


^iv{E-Vi){2T)-^ _ ^iv{E-Vi){2T)~^ 


E - Vi dv 

both times. Using {NT)~^ < 1, we bound the resulting integral by 


(5.19) 




G'^gs.i 

{V^ - U)2 Jo 


pOO ro 

I '■“I 


du 


1 


1 + \vr + |u 




Y\i\l + uci + iu(i/|V2 
\dj\ \dj 


UCi 


ivdi 


E 

j,k 


141 


1 + ucj + ivdj I 11 + uck + ivdk \ 


Define the random set of indices 

Es:= {I : \Xi{h!Y^) - E\ < N^-^. 

Note that by rigidity we have that 

E[|F 5 |“] < 

for any M > 1. 

By the Schwarz inequality we have 


(5.20) 

(5.21) 

(5.22) 


lg.M\Gn\^] ^ 


Gigs,. 
(Vi - U)2 


rCC nO 

I 


du 


Yli\l + uci + iud/|V2 
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1 + ucj + ivdj I 11 + uck + ivdk \ 


2 ^ 

j,ki$Ts 


14 


141 


+ ucj + ivdj I 11 + uck + ivdk \ 


(5.23) 
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The argument above leading to the bound (5.18) yields 


E[^i] < 


{Vi - EY 




(5.24) 


i.e., we bound Ai by the RHS of (5.16) except that we have a prefactor of (V) — E)~‘^ instead of T~^. 
The term A 2 can be bounded similarly; to be more precise, it is bounded by the same thing as (5.16) 
(with of course replaced by (V) — E)~‘^) except we take 4 additional d^/s for the denominator to 
deal with the in the numerator. We get 


E[T2] < 


1 


{Vi - EY 

We next handle A 4 . Define the random variable 

1 v-i 1 




Qs ■= 


N 


I \(d pi 


(5.25) 


(5.26) 


We have the following bound which is a consequence of rigidity and the proof of Lemma 7.5: 

E[Qf ] < (5 27) 

We provide the proof of this in Lemma 5.3 below. Repeating the argument as before we then get that 


E[2l4] < 


1 


EiQlAf] < N 


{V - EY 


10(5 


(5.28) 


where we have used the Schwarz inequality in the last step together with (5.17) and (5.27). 

The double sum in ^43 contains only \Es\‘^ terms. Since we have (5.22) we will just bound each 
term individually. We have, 


rCC rO 


du 


1 


\dj 


141 


— ^Gs,i 


WiY + uci + ivdi 1 1/2 11 + ucj + ivdj I 11 + uck + wdk\ 

p . 1 

Jo 


du 


£ 


du 


(1 + C„iU)V2(l + Ca2uY{l 

I dj dk I 

(1 + |4u|)V2(l + |4u|)V2(l + \di3v\) 


=: Ajk 


(5.29) 


Observe that due to the extra terms appearing in the denominator of the djdk terms of (5.23) it is 
possible to make the choice ai = j or k in the above inequality. This will be important later when 
we deal with the event {Nx^ > 2}. Let us first do the u integration. The region of integration is split 
into two parts: [0, 00 ) = [0, u 00 ). In the first region we bound the integrand by 1. The latter 
integral we bound by 


£ 


1 


--20r 


d/2. 


d/4 1/4 ,1+r 


dtt < 


42 „ 1/4 1/4' 


(5.30) 


Performing a similar split for the v integration we get the bound 


--20r 


Ajk < {s + 


d/2. 


..1/4 1/4 


)(^ 


10 


“'CKl ^02^j 


14 r 


) < e^N^Af + e 


-50r-l/2 a4 141 ^^^ 141^'^^ 

4 4 

We have used the definition of Qsj to choose Cq. > e/{2Af) and |4I ^ (2Ai)~^. Note also by the 
definition of Es that Cj,Ck > eN~^. We have used this bound for Cj and Ck, as well as the bound 
\di\ < s~^ which holds for any I, to bound the other terms that arise when the LHS of (5.31) is 
expanded (i.e., the terms of the form e40|4|i/2|4|i/2^|^^|?-^ etc.). Observe that 


\dl 

Cl 


< 


(5.32) 
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for any 1. We have therefore derived the bound 


E[A 4 ] < 


2(5^-l-50r 


{Vi - EY 


E[|F 5 |"(Af)] < 


8(5^-l-50r 


{Vi - EY 


where we have used the Schwarz inequality and (5.17), (5.22). 
Collecting everything we get 


mgsJGnl"] ^ 


1 


^115^-l-50r 


{Vi-EY + T^^ 

i.e., we have derived (5.15). As a consequence of (the proof of) Lemma 7.5 we have 


-S 

N ^ 


< nY 


N^ ^/{vY^eYTt^ 

Therefore after redefining 5 and r and using (5.14), (5.34) and (5.35) we get 


^ < Ce= (e^l^ + 1 f] 


j^S^-l/ 2 -r 




N ^,^iV,-EY + T\ 




and so we obtain the bound in the case k = 1. 

We now begin the k = 2 case. The starting point is the inequality 


(5.33) 


(5.34) 


(5.35) 


(5.36) 


P[iVx, > 2 ] < (Im [m^])^] < — ( 5 . 37 ) 

(i) 

where denotes the number of eigenvalues of the ith minor of Ht lying in the interval X^. We 
have used the interlacing property which implies that if two eigenvalues lie in then at least one 
eigenvalue of any minor lies in X^. 

Define the event Qs^i as before. Due to (5.14) we have 


On the event we will prove the inequality 

E[la,,l,„(.),i,|GiiP] < N\{Vt - E? + 


7 ’2 1 


]) 


l/(l+r) 


(5.38) 


(5.39) 


Proceeding in the k = 1 case we use the identity (5.8) and take the expectation over the ith row. Note 
{i) 

that is independent of the ith row of Ht- We see that 


1 r°° r 


q-^'^-^-'>^ei( 2T) ^+iv{E-Vi)(2T) ^ 
Y\iil + uci + iudz)V2 


(5.40) 


Again, we get two bounds on the integral and minimize over the two at the end. The arguments are 

(i) 

very similar to the k = 1 case, except that due to the fact that > 1, we will always be able to 
choose an index ai so that Cq,^ > . In the case A: = 1 we could only choose Cq,^ > e/(2A|) ~ eN~^ 

by rigidity. This was one of the key observations in the work [11]. 

For the first inequality, we proceed as in (5.16) and derive 


c r® 


0 |1 + CaiU|V2|l + + C, 


''«3 


u\ 


£ 


du 


(1 + |d/3iu|)V2(i + \dp^v\Yl‘^{l + \d0^v\y 


(5.41) 
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Since at least one eigenvalue of the minor lies in so we can choose As before, we then 

choose the other indices so that ^ £/( 2 A|) and |(i^J > ( 2 Ai)“^ which is possible by the definition 
of Gs,i- Hence, 

]E[lg, > l]V(i+-) (5.42) 


2 ^ 2 " 


where we have used the Holder inequality in the last step and (5.17). For the other bound we proceed 
as before and integrate (5.40) by parts twice in v. We obtain 


+ A 2 + A 3 + A 4 ) Bi + B 2 + B 3 + B 4 


(5.43) 


where the Aj are defined as in (5.23). The same argument leading to (5.42) (i.e., the same arguments 
we used to obtain (5.24) and (5.25) except taking Cq,^ > e~^) yields 


E[Bi + B 2 ] < 

Similarly we get 


C 


(Vi - EY 


e-’'E[A5l 


{7V«>i}J ^ (y._EY^ 




C 


jyl25]p|-^W ^ yl/(l+r)_ (- 5 ^ 44 ^ 


E[H 4 ] < 


< 


{Vi-EY^ 

^125^-r]p|-^0) ^ ;^jl/(l+r) 


{Vi - EY 


(5.45) 


using (5.27), (5.17) and the Holder inequality. We proceed as before in our estimation of H 3 . We 
estimate the general term by 


^Qs,i 


{NY>1} 


J 'OO rO 

VY- 


du 


1 


\dj 


WiY + uci + ivdi 11/2 11 + ucj + ivdj I 11 + ucfc + ivd^ 


*^^(1 + Cq,jU)V2(1 + Ca2lt)'’(l + «Cj)V4(l + 4(Cfc)H4 

^ r1^.41 

J — 0 ( 


(l + |d,u|)V2(l + 14^1)1/2(1 + |d^^|)r 


= : B 


(5.46) 


Note that we can take ctj = j or k due to the extra factors appearing in the denominator below the 
dj and dk terms in the first line of (5.46) (this was the observation made after (5.29)). Due to the 
fact that > 1 we can take Ca^ > ■ By the definition of Qs^i we can take Cqj > e/(2A^) and 

ha > (2Ai)“^. Splitting the two integrals into the regions [0, u [e^°,oo) as in (5.30) we get, in the 
same fashion as (5.31), the bound 


-- 20 r 


o, 1 /,io I _w no , 

jk — Gs,i {Af^®i>l}f^ 1/2 ^ 1/4 1/4^1^ IJ Ir ' 

= Cai Ca2^j 


\dpy 




--50r a4 


:5/V45^4 ^ g-50r+l/2^4 l^il^l4l^ ^ 


„V4 1/4 

S’ ^k 




<e 

From Holder and (5.17), (5.22) we get 


E[H3] < s-^^''N^^E[Af\Fs\'^l (i) ] < e-^°^/Vi25p(/Y0) > i]i/(i+r)^ 

i-^Xg —^ 

Collecting everything and optimizing over our two bounds we get 


(5.47) 


(5.48) 


^ “ (Hi - EY + 


N^e-^¥[N^"'> > l]V(i+D 


< 


1 


{Vi - EY + T2 


^41-3r 


(5.49) 
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where in the last step we have applied the k = 1 result to the ith minor of Ht- Plugging this and 
(5.38) into (5.37) and summing over i,j yields the claim. □ 

Outside the interval XEQ^qC we still have the following bound. It is a corollary of the above proof. 
It may also be established using the methods of [6,11]- 

Corollary 5.2. For all r > 0 there is a Cr > ^ so that 

¥[Nx, >k]< _ Eo\ + if- (5.50) 


for k = 1,2. 

Proof. The proof is essentially the same as the proof of Theorem 5.1. We illustrate the differences 
in the case k = 1. Since now the event may have probability 1, the bound (5.14) becomes 

E[lgfGuf<eN^. (5.51) 

On the event Gs,i there is no need to integrate by parts as we do not need to optimize the prefactor. 
We just use the argument leading to (5.18) and get 

mgs..\Gnf ^ (5.52) 

Since all the eigenvalues of Fkf lie in the interval [Av,N^^] + [—|11, |11] it is easy to see 
that 

E[Af] < N^^{\E - Eo\ + if^. (5.53) 

This gives the k = 1 bound. The k = 2 bound is similar. □ 

Proof of Theorem 3.6. We decompose the event < e/N} as follows. 

{|Ai,t - < e/N} = {|Ai,i - Xi+f < e/N, |A,,t - ^ 

00 

[J{\Xi,t - Xi+f < £/N,\Xi,t - 7 i,t\ >NfN,n-l< \\Ht\\ < n} (5.54) 

n^l 

= (5.55) 

n 

Define 

■= bi,t + e(j - 2)/iV, 7 i,i + e{j + 2)/iV]. (5.56) 

Using Theorem 5.1 we have 

P[7'i] < 2 ^ 2 ] < (5.57) 

\j\<N^/e 

Let Pi, p 2 and ps satisfy pf^ ~kP 3 ^ = 1 with pi > 1. We have by Holder’s inequality and Corollary 
5.2 


FiGn] < 2 ^ 2, \K,t - 7i,t| > NfN,n- 1 < ||IT|| < n] 

bj<4n/e 

< (P[||IP|| >n- e-^iog(iV)«/P3_ 


(5.58) 


For large enough N depending on p 2 and ps and Cr we have due to the exponential decay of the 
entries of W, 

00 

2 (P[||IT|| >n- 1 ]) Vp 2 g-i^iog(iV)e/P 3 < i ( 5 ^ 59 ) 

n^l 


Hence choosing pi close to 1 we obtain the theorem. 

We required the following lemma in the proof of Theorem 5.1. 


□ 
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Lemma 5.3. We have. 


E[Qf ] < ^5 gQ^ 

for E e lEo,qG o-iT-d Ts as in (5.21) and Qs as in (5.26). 

Proof. Let 1 > > g > 0 and define Qi to be the event that \^j^T — for j e .4.gj,T and 

||VL^*^|| < 3. Then Qi holds with (^, z^)-high probability. On the complement of Qi we just use the 
bound IQ^I < and so we may assume that Qi holds. We write 


N\ Fl 


= 1 


Si 


E 


1 


JV|Af - E\ 

■yj,TelEQ,q.^G 


LQi 


E 


1 


jt?, JV|Af - E\ 

lj,T<f^EQ,q^a 


= :A 2 +Ai. (5.61) 


For the terms in A 2 we get directly from the definitions of Qi and Ts that 


A 2 < 




|fc|<7V 


< N‘ 
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(5.62) 


Fix a q 2 satisfying qi > q 2 > <?• Define the intervals for k > 0 

=Eo + [-2'^Gq2, 2’^Gq2] (5.63) 

as in the proof of Lemma 7.5. By definition, on the event Qi the norm of is less than 3. We have 
the bound |Aj(T) — Xj(B)\ < HA — B|| which holds for matrices A and B. Hence we see that on the 
event Qs,i, 

\\j{H^^^)-\j{V)\<VT. (5.64) 

Therefore, for /c > 1 we see that 


|{i : A,(H«) e < |{i : H* e X«,}| < GN2^G 


(5.65) 


where the last inequality is proven in Lemma 7.5 and holds as long as k < |log(cG)/log(2)| for some 
c > 0. From this we see that 


■4 E 


CMC) (jj^2'^G 


fc=i 


2^0 


< Clog(iV) < 


The claim follows. 


(5.66) 

□ 


6 Proofs of main results 

Proof of Theorem 2.5. Fix t = N'^/N e %} and ko e Aq^t- Let 


to = t- 


N‘^ 


N 


( 6 . 1 ) 


for oj' < r/100. We apply Theorem 4.2 with the matching GOE matrix Htq = aW + b for a, b and jo 
chosen as in Lemma 4.1. It follows that there is an e > 0 and an event E satisfying PHT”] > 1 — N~^ 
on which 

N-^ 


^j+ko,t ^j'+ko,t (Aj+io,t k‘j'+jo,t) 


< 


N ’ 


( 6 . 2 ) 


for |j| + |j'| < for some c > 0. We take c < u/20. Recall that A are distributed as A(L7t) and p, 

are distributed as X{atW + bt). Clearly at satisfies |a — at\ < |t — to|- 
By Lemma 7.6, 


\Pic,to{'yko,to) Pic,t{Xko,t)\ ^ \Pfc,to{'yko,to) Pic,toilko,t))\ + |Pfc,to(7A;o,t) Pic,ti'yko,t)\ 

< log(iV)^^-^ < ^ 


to 




(6.3) 
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Since the eigenvalue differences appearing in (6.2) are bounded by N'^^/N « with very high 

probability, we see that there is an event J -2 with P[J^ 2 ] > 1 — on which 




at 


< 


N- 


N 


(6.4) 


for \j — ko\ + \j' — ko\ < N^, after possibly decreasing e. Now note that is distributed as the 

eigenvalues of a standard GOE matrix. The claim follows by expanding the test function O on the 
set j^ 2 - n 

Proof of Theorem 2.4. The above proof establishes the fact that there is a c > 0 so that for any 
index i e Aq^t and positive integers ki, ...,kn ^ we have 

|E^*[0(pfc,t(7M)(Ai - A i+fcl)) Aj+fc^))] 

-E‘^‘^®[0(psc(Aii)(Ai - Ai+fcJ, ...,psc(A*j)(Ai - Ai+fc„))]| < N~'^. (6.5) 


The method of going from the above estimate and optimal rigidity results to averaged bulk universality 
is standard and can be found for example in Section 7 of [13]. We omit the details and only mention 
that our result is stated with an upper bound on the size of the averaging window which ensures that 


- pi,,tiE')\ < N-'^ ( 6 . 6 ) 

for all energies satisfying \E' — E\ < b. If a larger averaging window is desired then p{c,t{E) must be 
replaced by pic,t{E') in the theorem statement. □ 


7 Proof of local deformed semicircle law 


7.1 Stability estimates 

In this section we will derive various qualitative properties of the deformed semicircle law. All state¬ 
ments and proofs here are deterministic. It will be notationally convenient to work with the definition 


and 


Ht = v + Vtw 


^ y TT- ^ -rw 

N f^^Vi - z - Tm{c{z) 


(7.1) 


(7.2) 


where we have dropped the subscript T in mfc,r- Throughout this section we will assume that V is 
(.^, G)-regular at the point Eq. 

We will need the following elementary lemma about Stieltjes transforms. 

Lemma 7.1. Let m be the Stieltjes transform of a probability measure p with support in [—A*, A*]. 
There is a eonstant C > 0 s.t. 


\rh{E -I- ir])\ < C(| log(? 7 )| -I- | log(A*)|) sup Im [m{E + ip')]. 

r]'>rj 

Proof. Define the dyadic intervals for j > 0 by 

Ij := [E - 2^p, E + 2^p] 

and set X_i = 0. We have 

2 ^p 


Im 


[rh{E + 2^p)] > J 


r {x - EY + (2iry)2 


dp{x) > c 


2^p 


(7.3) 


(7.4) 


(7.5) 


Let = supAj : p{Ij\Ij-i) ^ 0}. Then < C{\ log(r/)| -h | log(A=,)|). We have 


n 

\m{E -I- ir/)] < V _ 


Mx) ^ ^ 




dp{x) 


< 




p(Ij) 


'^ < Cj^ sup \m{E + ip')\. {7.6) 


where we applied (7.3) in the last inequality. The claim follows. 


□ 
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Lemma 7.2. Assume that V is {i,G)-regular at Eq. For 1 > q > 0 and oj > 0 all of the following 
hold for E £ lEQ,qG o,nd 10 > r; > N~^ and T e %], for N large enough. Below the constants do not 
depend on the choice of uo or q. For the Stieltjes transform, 



c < Im \midzy\ ^ C. 

(7.7) 

and therefore 

cT <\Vi — z — Tmiddl ^ C*. 

(7.8) 

For the density we have 

c<pidE)<C, \ddE)\<^. 

(7.9) 

Finally, 

c<\1-TR2{z)\<C, \T‘^R3 {z)\<C 

(7.10) 

where 

Rk{z). ^ y ^ ' 

(7.11) 


Proof. By the assumptions on V and Lemma 7.1 we see that there is a constant Cy > 1 s.t. 


lmv(E + i? 7 )| < Cv log(iV), E e Zeo,g, V > T/{lQCv) (7.12) 

and 

Cy^ < Im [mv{E + ir/)] < Cv,E e 1 eq,g 10 > ry > T/(10C'y). (7.13) 

For notational simplicty we dehne the set 

A. = {E + ir] : E £ Ieq,Gi W > rj > T/(10Cv")} . (7.14) 

Fix e 'lEo,qG and define 

rji,: = inf {r/ < 5 : \mic{E + ir?)| < 2Cv log{N), and {2Cv)~^ < Im \mic{E + irj)] < 2Cv} ■ (7.15) 

Using the bound |mfc| < 1 for r/ > 1 we see that rj^ < 1. For a contradiction suppose that rj^ > 0. 

Since mfc is continuous in the upper half plane and the set defining ry* in (7.15) is nonempty (it has 
at least all ry >) one of the following three possibilities must hold. 

(i) |mfc(F^ + iry*)| = 2Cv log(A^) and {2Cv)~^ < Im [m{c{E + iry,^)] < 2Cy 

(ii) |mfc(F^ + iry*)| < 2Cy log(A^) and (2Cy)“^ = Im [mfdE + iry,^)] 

(hi) |mfc(F' + iry*)! < 2Cv log(A^) and 2Cv = Im [mfc(S + iry*)]. 

We have the equation 

mic{E + iry^f) = mv{E + ir]^ + Tmic{E + iry*)). (7-16) 

Suppose that (i) holds. For N large enough depending only on Cy, we have that i? + i?y* + Tm^dE + 
iry*) e A. But then by (7.16) and (7.12) we have 

|mfc(F; + i7y=,)| < C'ylog(7V) (7.17) 


which is a contradiction. If (ii) holds then still + iry,^ + Tm^dE + ir/d £ A. By (7.16) and (7.13) we 
get 

Im [mfc(F; + iry*)] > {Cy)~^ (7.18) 

which is a contradiction. If (iii) holds we arrive at a similar contradiction. Therefore, ?y* = 0. We 
conclude (7.7), (7.8) and the Hrst bound of (7.9). 

Taking imaginary parts on both sides of (7.2) and rearranging we obtain 


/ T ^ 1 

Im rmfe(z)l 1-y ^- 

L tCV ^ rr _ . _ 'I 




= h 


1 ^ 

-y 

N Aj 


N ^.^\Vi-z-Tmidz)\^ 


(7.19) 
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and so 


N 


1--V-- ^ 

N ^ V, - 


Denote z + Tmic{z) = a + 5i for a, 6 e M. We have, 

rr. N 


Re[l - TR 2 {z)] = 1-^2 


1 


> 0. 


rr N 

2b^-y 

N ^ 


N^,\V,-a-bi\^ N^^\Vi-a-bi\^- 


(7.20) 


(7.21) 


By the hypothesis on V and the fact that a e Teo,g we see that there is a C > 0 so that if 6 > CT 
then 


T ^ 1 

iV^ \Vi-a-bi\^ ~ 2' 
2=1 ' ' 


(7.22) 


We have by Lemma 7.4 below that there are at least NT entries of V less than distance CT from a. 
Hence in the regime CT > b > cT we have 


. T ^ 1 NT 

— > - > -> c 

~ NiCT)^- 


(7.23) 


and we conclude the first part of (7.10), using (7.20) (note that by (7.7) we always b > cT). The 
second part of (7.10) is an obvious consequence of (2.3). The last thing to prove is the second bound 
of (7.9). Differentiating (7.2) gives 


N 


dzmc{z) = 1-^2 


-1 N 


^ {Vi-z- Tm,,{z)f j N {Vi-z- Tmfc(z))2 ' 
Therefore, by the first bound of (7.10), 


(7.24) 


N 


\dzmic{z)\ < C— y 


N \Vi- z- Tmfc{z)\ 


C 

2 — 'IP 


(7.25) 


where in the second inequality we have used (2.3) (note that we have used the fact that |Re[mfc]| < 
T~^/^ which implies that if Li e ^Eo,qG then E + TRe[mfc(Li + irj)] e Teo,-^g)- We conclude the bound 
for by the Stieltjes inversion formula. □ 

Lemma 7.3. For z e T >2 where T >2 is defined in (3.11) we have, 

\Vi-z-Tmifiz)\>c, c<\1-TR2{z)\<C, \T'^Rz{z)\<C (7.26) 


and 

\mifiz)\<C (7.27) 

Proof. The first bound follows from the fact that Im [^] > 10 on T> 2 - The other two bounds are 
consequences of the fact that we have \Rk\ < C for z e 'D 2 - The bound (7.27) is obvious. □ 

For the above proofs we required the following elementary result. 

Lemma 7.4. Fix 1 > g > 0 and uj > 0, and V be {i, C)-regular at Eq. Then for N large enough we 
have uniformly for E e lEo,qG oind p ^ Tu, 

cNrj < |{i :Vie[E-ri,E + rj]}\ < CNrj (7.28) 

for some eonstants depending only on the constants in (2.3). 
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Proof. Let = \^E — r],E + rj]. The upper bound follows immediately from the fact that 

Vi sir, {Vi - Ef + r/2 < 2if. (7.29) 

Denote by C* the constant appearing in the upper bound just obtained for (7.28). For A: e Z denote 

= + (7.30) 

and define 

(7.31) 


ko = inf |A:| : n {R\Ie,,^,g} ^ 0 


We have, 


1 


Im [mv{E + irj)] < S V S 


V 


Tf 

^ {Vi-Ef + 

\k\<ko 


(7.32) 


Note that C'qjG'^ < CN ^ by assumption on rj. Fix X > 0, L G N and let Nl := |{f : V) e 
[E — Lt], E + Lri]}\. We have, 


L V L 


V 


< 


L V L 


r] 




V 


N E)^ + ~ N ^„(Vi-EV + ri‘^ N ^,^(Vi-EV + T]'^ 

\k\<ko ViEli,^ kQ>\k\>L ViSZ^^ 


^ Nl ^C* ^ Nl 20C* 

- L L ■ 

I /c I !7> X/ 

Choosing L > 100C'*/c\/ where cy is the lower bound from (2.3) we get for large N, 


Nl > (Nv)- 


which yields the claim. 
Define 


gi{z) = 


Vi- z- Trufciz)' 

We will require the following bound later. 

Lemma 7.5. For all z e EL^q, we have 

< ClogN. 


(7.33) 

(7.34) 

□ 

(7.35) 


(7.36) 


Proof. This is clear for z e 1)2 ■ For z e Vi we write 

L i 9 *wi + j^ E tewi=:A 2 + y. 

I i-.VielE^yqQ i-.Vi^lEQ^yqC 

We first bound 2 I 1 . Define for k > 0 

Xi') := Eo + [-2"v^G, 2^VgG], :=\{i-.Ve 


(7.37) 


(7.38) 


For k < ki := [—log(G)/log(2) + Ij we have by the assumptions (2.3) and the proof of Lemma 7.4 
the bound 

N^^'^ < C2^GN. (7.39) 

Hence, 

^ > > \ n.-17.) \ -L 

N 


l<fc<fci y,gx(l)\x{l)^ 




< 




N, 


( 1 ) 


G 


N ^ 2>^G 

l<k<ki 


G 


<Gki + G <Glog{N). 


(7.40) 
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For the term A 2 we define the intervals for k e Z 

:=E + rRe[mfc(z)] + 2Tk + [-T, T] 

and let 

wf i-KiiV'.elfll. 

( 2 ) 

Let k 2 and k^ be the index of the leftmost and rightmost intervals If, ^ that intersect lEo,y/qG, 
respectively. By Lemma 7.4, 

< CTN, for k 2 <k<kz. (7.43) 

Then, 


(7.41) 


(7.42) 




k2<k<k3 y.gi(2) 


< 




N, 


( 2 ) 


C 




< ClogA^. 


(7.44) 

□ 


We also require some results about the qualitative properties of the time evolution of mfc,t. 
Lemma 7.6. Let mfc,t satisfy (3.5). Then 


We also have 


dtmic,t{z) = ^d^[mfc,t{z){mic,t{z) + z)]. 


= -Re[mfc,i(7i,i)] - 


Fix cj > 0 and q > Q. We have uniformly for i e Ag^to CLnd t,to ^ Tuj, 

< ClogiV. 

as well as 

\dtPic,t{E)\ < j 


(7.45) 

(7.46) 

(7.47) 

(7.48) 


for E e lEo,qG- 

Proof. Equation (7.45) is a straight-forward exercise differentiating (3.5). The equation (7.46) follows 
as in the proof of Lemma 4.4 of [22]. The resulting bound (7.47) follows from the bound [rufc^tj < 
ClogA^ which is obtained from Lemma 7.5. The bound (7.48) follows from (7.45) and the fact that 
we bounded \dzm{c,t\ ^ C/t in the proof of Lemma 7.2. □ 

Proof of Lemma 3.4. This follows directly from (7.47). □ 


7.2 Local deformed semicircle law 
7.2.1 Preliminaries 

We collect here some identities and tools required for our proofs. For a set of indices T [1,1V] we 

(T) 

denote the minor of Ht obtained by removing the columns and rows in T by Hf '. Similarly we denote 

fT) (T) 

the Green’s function of Hf ^ by GL '. For the minors we will still normalize the empirical Stieltjes 
transform as 

= ^tr I ^ (7.49) 

li rp Z 

We use the notation Tf = T u {f} and for T = {i}. We use the notation 

(T) 

S - S (7-50) 

ij 

m 
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to denote summation over all indices not in T. 

For a row i and random variable X we denote and = X — Pj[X]. 

In the following lemma we collect the resolvent identities required for the proof. We refer the 
reader to, for example, [9] for a proof. 


Lemma 7.7. The Schur complement formula is, for i ^ T.' 






For i,j, k ^ T and i,j ^ k we have 


^(T) ^ ^(Tfc) , '^ik '^kj 1 


^ik '-'ki 




G 


(T) 

kk 


UT) UTk) ^(T)^(Tfc)^(T) 

'^kk 


For i,j^T and i j, 


and 


(Ti) 


(Tj) 

^(T) _ _^(T) Y* h - -r’W V 

'^ij ~ '^ii 2-1 ~ ^jj 2-1 '^ik '^kj 

k k 


Gij) 


qG) _qG)qG’'')( f^.. _ y* /j. .) 

^ij ^ii ^jj vG 2-i "'^‘^^rnn "'nj)- 

m,n 

We have also the Ward identity 




(T) 

a * 


(7.51) 


(7.52) 

(7.53) 

(7.54) 

(7.55) 


We collect here the large deviations estimates we will require. We refer the reader to, for example, 
[18] for a proof. 

Lemma 7.8. Let (a*) and [hi] he centered independent random variables with variance cr^ and having 
subexponential decay 


P(|ai| > xa) < P(|6i| > xa) < 


(7.56) 


for some positive constant Cq and 6 > 1. Let y e C and Bij e C. Then there exist constants oq > 1, 
Aq > 10 and C > 1 depending on 6 and Cq s.t. for oq < ^ < log log and ip = {logN)^, 


N 


Xj 

J 

O' Bi> 

i=l 
N 

P I ^ diBijOj 


N 


i=l 


/ N \ 


< e 


-(logAf)« 


TV \ V2^ 

|2 
\ 


> if^o ( ^ \Bii\'^ j I < e 


-(logAr)« 


',*=1 




N 

X diBijbj 
*J=i 






TV \ V2^ 

2 

•ij\ 

Wj / 

TV \ V2N 

2 

■ij\ 

2=1 


< e 


-(logAr)« 


< e 


-(logAr)« 


for N sufficiently large. 


(7.57) 

(7.58) 

(7.59) 

(7.60) 
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7.2.2 Proof of Theorem 3.3 

The proof follows closely the proof of Theorem 2.10 in [21]. We will just note the differences here. 
The first difference is that we consider V fixed instead of random. This change was already dealt with 
in the proof of Theorem 3.3 of [22]. The extra randomness in V corresponds to the larger error terms 
which appear throughout [21]. For us, taking V fixed corresponds to setting A = 0 everywhere in the 
estimates appearing in Sections 3 and 4 of [21]. Additionally, in the paper [21] the local law is also 
proved at the edge of the spectrum. We will not deal with the edge case and will therefore avoid some 
of the technicalities encountered in [21]. We split the proof of the local law into two parts, first for 
z e Vi and then for z e 7^2 ■ The proof in the domain Di is similar to [21] and we now summarize the 
changes. 


7.2.3 Local law in the spectral domain Vi 

The major difference between the proof here and that in [21] is of course that we consider Ht = 
V + ^/TW for small T whereas in [21], T is order 1. We summarize what difficulties this causes and 
how they are dealt with. The first lies in the stability estimate (7.8) which only gives a lower bound 
of cT whereas the analogue in [21] is Corollary A.2 where the lower bound is c. This changes a few 
steps of the proof. Instead of (3.17) of [21], the Schur complement formula is written as 


Gii = 


where 

and 


Vi - z - Trufc + T'tp + Yi 
'll; = rriic- tun, Yi = VTuuu + T{mN - - Zi 


Z, = Y,hrkG^Shu-Tm^S = Qi 


(d 


kl 




k,l 


(d 

with 'nvJ the empirical Stieltjes transform of the ith minor of Ht- 


Let 


Define 


A{z) := I'll;] = \mN{z) - 'mic{z)\, 7f{z) := 


9i{z) ■■ = 


Vi - z- Tmfc{z)' 

We have the following a priori bounds on Zi and Yi on the event that A( 2 ) 
version of Lemma 3.8 of [21]. 

Lemma 7.9. O'n the event {A( 2 ;) < 'we have 'with {^,u)-high probabilit'y, 


(7.61) 

(7.62) 

(7.63) 


(7.64) 

(7.65) 

o(l). This is our 


maxjZj] < CT7>, maxjgjYi] < 


(7.66) 


and 


max I Kj < C 


I Vn 




(7.67) 


for a constant C. 

Proof. By the large deviations bounds (7.58), (7.59) we have with (^,zz)-high probability 


k,l 


1/2 


< GT7/ 


(7.68) 
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where we have used (7.55) and the inequality [21] 

\mN - < {Nr])~^. 


(7.69) 


The bound (7.67) for Yi follows from this and the subexponential decay assumption for wu which 
yields with (^, i/)-high probability, 

«.«l Vf^. (7.70) 

For the bound (7.66) for Yi we use (7.69), (7.70) and \gi\ < min{r/~^,T“^}. □ 

Using these estimates, we can prove the following analog of the weak-self consistent equation which 
is Lemma 3.9 of [21]. 

Lemma 7.10. On the event {A(z) < we have with {^,v)-high probability, 



(1 - TR2{z))i; - T^R3{z)'iP^\ < CT + C:^ 


(7.71) 


Proof. By Lemma 7.9 we have on the event {A( 2 ;) < ip that \Yi\ < To(l). Therefore, by the 
stability bounds (7.8) and (7.26) we can expand (7.61) and obtain 





V L + V L 9? +r,) + L L 

i i i 

1 2 afmYi + Y^) + O ( (max + T^A^))N-^ ^ \9J 

i \ j 


(7.72) 


We use a simple Schwarz inequality < C'|yi|^(log A^) -l-log(A^) ^A^ on the first term in the second 
line of (7.72). By Lemmas 7.9 and 7.5 we have 

<C'log(iV)T^ (7.73) 

i 

for r > 0. This, together with (7.66) yields the claim. □ 

With this, one can derive the weak local deformed semicircle law. 

Theorem 7.11 (Weak local deformed semicircle law). The event 


n 

TsTlo 


|mAr(z) - mfc(z)| < C 






(7.74) 


holds with {^,h')-high probability. 

This is derived from the weak self-consistent equation exactly as Theorem 3.1 is derived in [21]. 
To be more precise, one proves the initial estimates for g = 2 va the same fashion as Lemma 3.10 
of [21]. The dichotomy argument of Lemma 3.12 is easier as we have that the LHS of (3.45) of [21] 
can be taken to be 1, and so one only needs to do the ‘Bulk case’ argument in the proof of Lemma 
3.12 of [21]; i.e., we will only be considering the regime g > g(U,E). The stability of the coefficients 
in the self-consistent equation is provided by the estimates of Lemmas 7.2 and 7.3 (as we are in the 
bulk the coefficient (1 — TR 2 {z)) does not degenerate which means that a > c where a is defined as 
in [21]). 

The proof of the weak local deformed semicircle law is then completed as in [21]. We record here 
the following corollary of the above proof. We will not need it in this paper but record it for possible 
future use. 


Corollary 7.12. Suppose that one has that c < Im [myiE + ir;)] < C for g > go where go satisfies 
(for example) N~^ < go < N~'^, and all E e T where I is some fixed interval. Let gj = p^^^^go for 
j = 1,2. Then with {f^,v)-high probability, one has 


\mN{z) 


mfc,T(2:)| < C 




(7.75) 


for all 0 < T < N 


, E eX' and g 2 < g < W where X' is an interval compactly contained in X. 
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Proof. It is easy to see that the above proof directly implies the weak law for T > r/i with T « 
for rj down to the optimal scale r] > N~^. On the other hand, the above proof also yields the weak 
law for 0 < T < r/i for 77 > 772 - For this, the key input is that the bounds (7.7), (7.8), and (7.10) 
hold in the regime rj > 772 . The first two are obvious as we have Im [rufc] = Im [mv{z + Tnifc)] and 
Im [z + Tmfc] > r ]2 » r/o- The upper bounds of (7.10) are clear as 772 » T. The lower bound is slightly 
nontrivial but we have, 


Re[l — TR 2 {z)] >1-Im [mv{z + TmfJz))] > - (7.76) 

V2 2 

for large enough N. The expansion in the proof of Lemma 7.10 is still possible as now one has the 
stability bound \vi — z — Tm{c{z)\ > 772 and iTip] + \Yi\ < To(l) + (T/A^)F2(^? < 7720 ( 1 ). 

□ 

We now return to the proof of Theorem 3.3. In order to get the strong self-consistent equation. 
Lemma 7.14 below, we require the following fluctuation averaging lemma. As the proof is somewhat 
lengthy, we defer it to the next subsection. We introduce deterministic control parameters 7 ( 2 ;) and 
<I>(z) by requiring 

7 ( 2 ;) < 93 “^^, 4 >( 2 ;) = < 73 ^^ ^ 

Lemma 7.13. Suppose that the event 

A{z)<j{z), \/zeVi (7.78) 

holds with {^,v)-high probability. Then, 

(^) 

with (^ — 2 , u)-high probability. 

The above lemma allows us to deduce the strong self-consistent equation. It is our version of 
Lemma 4.5 of [21]. 

Lemma 7.14. Suppose that the event 


< C'<y3i°«4>2. 


(7.79) 


A( 2 :) < 7 ( 2 ;), V 2 ; G Vi 


(7.80) 


holds with {^,v)-high probability. Then, 


A 2 


1(1 - TR2)^I; - < C- -- + 

logiV 

with (^ — 2 , v)-high probability. 

Proof. We proceed as in the proof of Lemma 7.10 and expand (7.61). We get 

V L c.. - ^ E »< + V L 9? + V L 


Note that we have 


Qi 




= Qi 


(i) 


VTwu - z-'^ hikG^^lh^ 


kl 


= Vfwii - Zi. 


(7.81) 


1 + G (^(max hfmf + T^A^))N-^ ^ ' 


(7.82) 


(7.83) 
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Hence, 


Yi = T{mN - + Qi 

Therefore, by Lemma 7.13 and (7.69) we have 








(7.84) 


(7.85) 


with (^ — 2, i/)-high probability. The other terms on (7.82) are bounded as in the proof of Lemma 7.10 
using the a-priori bounds in Lemma 7.9. □ 

The proof of Theorem 3.3 now follows from the strong self-consistent equation as it does in the 
proof of Theorem 2.10 at the end of Section 4 in [21]. 

7.2.4 Proof of Lemma 7.13 

Lemma 7.13 follows directly from the following moment bound and Chebyshev’s inequality. 

Lemma 7.15. If the event 


A(2;) < 'y{z), 'iz e Vi 


holds with {f^,v)-high probability, then 


E 


2r 


< (Cr)^^(<y93«$) 


4r 


(7.86) 


(7.87) 


for every r < i/(log A^)^“^/^/2. 

For its proof we will require the following a-priori bounds on the Green’s function matrix elements. 
Lemma 7.16. Suppose that the event 

K{z)<-i{z), VzePi (7.88) 

holds with {f^,u)-high probability. We then have uniformly for |T| < (logA^)^“^, 


9i 


a 


(T) 


- 1 


and for 0 < s < 1 and i ¥= j 


and 


\G\ 


m. 






25$ 




\9i 


Q 


* 

L 




(7.89) 

(7.90) 

(7.91) 


Proof. From the Schur complement formula we have 


9i 


G, 


(T) 


- 1 


= \9i 


Tf; + Y 


(T) 


(7.92) 


where 


(Ti) 


yG^ := Vfwn + nruN - ) - zP , := ^ h^kGG^ hu - . (7.93) 


kl 
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The proof of Lemma 7.9 extends to and uniformly for |T| < (logA^)^ ^ and we conclude 
(7.89). Moreover, we see that this observation also proves (7.91) after using 


Qi 






(7.94) 


For (7.90) we start with (7.54) and obtain with (f, z^)-high probability by the large deviations 
bounds, (7.89) and (7.55) 


ir’Wi _ ir^(T)^(Ti)| 


(Tij) 

h- — V h- 

/ j '^lYn^rnn '■ 
m,n 


nj 


^ ^^\9i9j\ 


Vn 


/im m 


(Tij) 

N 


Nrj 




In the last line we have used \gi\ < T ^ oi \gi\ < t] ^ as appropriate. 

Proof of Lemma 7.15. We illustrate the method by doing r = 1. We have. 


E 




G 


2=1 

= + A 2 . 

For the first sum we apply (7.91). Hence, 


ij 

1 ^ 1 

E Ift [0.7] I' + m L [7<3. [0.7] s^Qi 


07 


N 


Hi < 2 \9i? < ^ 


2=1 


where we have used 


N 


2=1 


*i'< 


We apply the second identity of (7.52) to the summand in A 2 . We have 


^[9^Q^ [G-y]9]Qj 


G 


-1 

33 


= E 


= E 


91: Qi 
9i Qi 


1 GijGji 

/^U) /=< /=v 0 )^ 

'^ii'^ii ^33 


GijGji 


9j Qj 


GuG, 


3^^^3 


/^(O /~fG) 

L^jj 


/^U) ^ 


9]Qj 

GjiGij 


1^33^jj \ \ 


where we have used 

E[XQi[Y]] = E[Q,[XY]] = 0 

for X measureable wrt to the ith minor of Ht- From Lemma 7.16 we see that 


(7.95) 

□ 


(7.96) 


(7.97) 


(7.98) 


9l 

< 9 

r’(Ti) 

'-'Ik 

rC^i) 


r’(T2) 

'-'kk 




G 


(TTi) 

Ik 


G 


(^ 2 ) 




I 7^ k 


(7.99) 

(7.100) 

(7.101) 


holds with (^ — 2, i/)-high probability uniformly for |Ti|, IT 2 I < (log A^)^ From these inequalities we 
deduce 


E 


glQ^ [Gu^ 9 ]Qj [ g -/]]| < \9i9j\{^^^^) 


(7.102) 
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and so 


(7.103) 


A 2 < ^ 

hj 

where we have used the fact from Lemma 7.5 that 

i 


This proves the case r = 1. 

For the general case we will just explain how to adapt the proof of (4.5) of [21] to our set-up. We 
start with (4.28) of [21]: 


E 


N 




1=1 


2r 


JV» 


E E l(r-r(i)|E[s?.(3„(G-;,.)„.4(3„„(G-;,„)] (7.105) 


rG7^2r 


where V 2 r denotes the set of partitions on 2r letters and r(i) denotes the element of 'P 2 r defined by 
the equivalence relation a ~ 6 iff Zq = and i = (ii, ...,Z 2 r). 

Fix now a partition F and let i satisfy r(z) = F. We apply the same algorithm as in [21] to the 
summand on the RHS of (7.105). We omit the precise details of the algorithm. The result is the 
expansion into the sum 


ngiQn{G-\0.-.giQr,AG-X)] 


(Ji ,...,f72r 


(7.106) 


where the ak run over the finite binary sequences generated by the algorithm and are the 

associated monomials in the resolvent entries. We will eventually apply (7.101) to the monomials 
{Fi^)„k ™ the above expression. To see that we can do so, we require the following lemma. It is an 
easy consequence of the definition of the algorithm generating the monomials which uses the identities 
(7.52). We give the full details of the proof. We encourage the reader unfamiliar with the algorithm 
to consult [21] or Appendix B of [8]. 

Lemma 7.17. Let {Fik)ak ® monomial generated by the algorithm outlined after (4.29) in [21]. Let 
F\ be the set of off-diagonal resolvent entries in the numerator and F 2 be the set of diagonal resolvent 
entries in the denominator. Suppose that Fi ¥= 0. Note that IJ- 2 I = +1- There exists an injective 

function tt : Fi —>■ F 2 such that, 

1. For e F\, the lower index of 7r{G^JJ) is either x or y. 

fT") 

2 . The unique diagonal resolvent entry Gxx £ .T2\7r(T'i) not lying in the image o/vr satisfies x = i^. 
That is, its lower index is the same as the index Qi^ over which the partial expectation is taken. 

Proof. Recall that the off-diagonal resolvent entries are generated by applying either of the rules 
in (7.52) to an existing monomial. Let ak = 0102 ...a„ for Oj e {0,1}. For each step j we construct 
TTj : J-'ij- ^ 02 ,j satisfying the lemma, where Fij and F 2 ,j are the off-diagonal and diagonal resolvent 
entries of {Fik)ai...aj, respectively. 

If the first step ai = 0 then there is nothing to do as there are no off-diagonal resolvent entries 
yet. If ai = 1 then 

= -(7.107) 

for some index j. We can take 

>ri(G.„)-G.«L >ri(G„.J=Oj,. (7.108) 

We proceed by induction. Assume that we have constructed -kj. First consider the case aj+i = 0. 
If the second rule of (7.52) was applied then we can let t^j+i = ttj if the rule was applied to a 
diagonal resolvent entry not lying in the image of vTj. If the rule substituted G^^'^ with 
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then we only have to define Tij+i{G^Jy^) = For all other off-diagonal entries 

we just let TTj+i = TTj, because the other elements of coincide with the other elements of J'lj 

and the other diagonal resolvent entries remain unchanged. 

(T) (Tl) 

On the other hand suppose that the first rule was applied to Gxy which became Gxy ^. We can let 
('TI'\ fT') 

Trj+i{Gxy ) = T^j{Gxy). The other elements of Tij+i coincide with the other elements of Tij- and we 
set vTj+i = TTj on these elements. 

Now we can suppose that Oj+i = 1. Suppose that the first rule was applied and that we have the 
replacement 

^(TO^(TZ) 

Gg) - (T)" ■ (7.109) 

^ii 

(T) 

Suppose that the lower index of TTj{Gxy ) was x. Then we set, 

’',+i(GT)-v(Gg>), (7.110) 

and set = ttj on the other elements of .Tij+i which coincide with the other elements of The 
case when the lower index of TTj{Gxy ) is y is similar. 

If the second rule was applied then we have the replacement 


G 


1 

(fy 

XX 


G 


(T)^(T) 
xy ^yx 


^xx ^xx ^yy 


(7.111) 


We then set 

7r,+i(Gg)) = GZy\ vr.+iCGg) = Gg> (7.112) 

and set t^j+i = iTj on the other elements of This yields the lemma. □ 

As a corollary we have 

(7.113) 

where h{ak) is the number of ones in the string a^. If b{ak) = 0 this follows from (7.91). If h{ak) > 0 
then note that the number of off-diagonal entries in the numerator equals b{(Tk) + 1. The bound (7.113) 
then follows from Lemma 7.17 and (7.101). 

For a label a e {l,...,2r} we let [o] denote the block of a in F. We let 5(F) := {a : |[o]| = 1} 
denote the set of single labels in F and denote by s = |5(F)| its cardinality. We wish to prove that 
any nonzero term on the RHS of (7.106) satisfies 




(7.114) 


First, suppose that one of the monomials, say is not maximally expanded (recall that a 

monomial is maximally expanded if each resolvent Gxy in the monomial satisfies 5(F) Q {a:,y,T}). 
By the definition of the algorithm generating the monomials, we must have that (F)j,)o-j, contains at 
least 4r off-diagonal resolvent entries in the numerator. Then h{ak) > 4r — 1 and so from (7.113) we 
have 


\9^,Q^AF^,)a,.■.9^,Mi2AF^,X,^ < (<^ 2 €^) 1+1 + .... + ( 4 .) + 1 ... < 

So we may assume that each monomial in (7.114) is maximally expanded. As in [21] we observe that 
for every single label a e 5(F) there is a label b e {1,..., 2r}\{a} s.t. the monomial {Fi^)„^ contains an 
off-diagonal resolvent entry with ia as a lower index. 

Hence, 

2r 

Yi K^k) > s (7.116) 

fc=l 

and we get the claim from (7.113). 
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For our fixed partition F denote its size by / = |F|. It follows from (7.114) that 
1 _ 1 N N 




2 l,... ,Zr 


(7.117) 

where di is the size of the iih block of the partition F. Using < 'q~^ we then bound the sum 


fci=l = l 
.-1 


iV 2 ^ 


TV TV , 

nr S -S U 


ki = l ki = l 


< 


Nr] 

1 

Nr] 


2r-l 2 ^ ^ 

E-E l^fci I 


N^ 


2r-l 


fcl=l fci=l 




log(A^)‘ < V 2 ^(^) 


(7.118) 


We have the inequality 2r + s + (4r — 21) > 4r. Inserting the two bounds just derived into the RHS 
of (7.105) we get 


E 


N 


V 2 »•'«■( TV 


i=l 


1 




2r 


< 2 ((^3S^)4r < (Cr)^’'(v93?^^ 

rGp 2 r 


4r 


(7.119) 

□ 


7.2.5 Proof of local law in P 2 

The proof of the local law in T >2 is similar but shorter than the proof in Di. The starting point is once 
again the Schur complement formula which yields 


mjy - rufc = 


1 ^ 

vE 


1 


1 


N Vi - z - Trriiciz) + T'ip + Yi Vi - z - Tmfc(z) 


2=1 

N 


-T^ - TYi 


ly 

N t[ {V-z- Tmi,{z) +T^ + Yi){Vi - z - Tm^z)) 


(7.120) 


From the fact that r/ > 10 on 'D 2 and that |yi| < with (^,z/)-high probability we obtain 

, , ^ \mN - mfcl , £ 1 

\mN-mfc\< - - - + (p‘^-^ (7.121) 

with (^, z^)-high probability and so Imw — mfc] < (p^f-^/N with z/)-high probability, uniformly in 'D 2 . 

We now proceed as in the proof of Lemma 7.14 and expand the Schur complement formula and 
obtain 


rriN = 


N 


1 

N 


TV 

^9i + 
2 = 1 
N 


'L9i + 

2 = 1 


L|;<,?(r^ + Tr,) + o('l|;i9,iv4] 

T =1 \ T =1 / 

1 E snTi> + ft[o.y]) + o (4 E isiivi]. 

T=1 \ i=l / 


(7.122) 


In the first line we used that |?/>| + \Yi\ < p^/\fN with (^,i/)-high probability, and in the second line 
that Yi = (5j[G4^] + 0{{Nr])~^). Since |1 — TR 2 I ^ c uniformly in V 2 this yields with (^, z^)-high 
probability 


Wn - m-fcl < 


g{z)ip^ 


N 


C 


1 ^ 

2=1 


(7.123) 


where we defined 


9{z) 


N {Vi - E)"^ + 7?2 


—Im [mv{z)]. 

T] 


(7.124) 
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Since g{z) < r] ^ we see it suffices to show that 

N 




N ‘ , 

2 = 1 


< 


Nrf 


(7.125) 


with (^, z^)-high probability uniformly in T> 2 . Before proving (7.125) we note that the proof of Lemma 
7.16 yields 


Lemma 7.18. Uniformly for |T| < (log A^)^ ^ and z eT >2 we have 


and for i ¥= j 


and 


-^-1 

“ Viv’ 

22 


r’(T) 


^ \9i9j 



Qi 




< 


^/n' 


(7.126) 

(7.127) 

(7.128) 


Proof of (7.125). We proceed as in the proof of Lemma 7.13. We use the same notation appearing 
there and will not redefine it. As in the proof of Lemma 7.13, we start by writing 


E 


N 




2=1 


1 


2r 


1 


2 S 

Y^'p2r ^1 )••• 


(7.129) 


We apply the same expansion algorithm as in the proof of Lemma 7.13 to each term in the above sum 
and estimate the resulting terms. With b{ak) denoting the number of ones in the string Ck, we claim 
that 

When b{ak) = 0 this follows from (7.128). For the moment let i{ak) denote the number of off-diagonal 
resolvent entries generated by the algorithm in the numerator. If b{ak) > 0, note that i{crif) > 6 ((Tfc)-|-l. 
Note also that there are i{ak) + 1 diagonal resolvent entries in the denominator. We want to use the 
bound (7.127) and the bound < 2\gi\~^ to conclude 

(7.131) 

The problem is that \gi\ can be quite small and so our estimate \Gii\~^ < 2\gi\~^ can be quite big, 
and we need to compensate for this with the factor \gigj\ that appears in the estimate (7.127) for the 
off-diagonal entries. 

For the monomial (^ 4 ) 0 -^ define Pj{{Fiif)aif) to be the number of times a diagonal resolvent entry 
with lower index j appears in Define to be the number of times that j appears as 

a lower index of any off-diagonal resolvent entry. For example, for the monomial 


^kk 


(7.132) 


we have s* = 2, = 2, p, = 2 and = 1. For (7.131) to hold, we need pt < Si for every i and every 

monomial. It is easy to see that by definition of the expansion algorithm that we will always have 
Pi < Si at each stage of the algorithm as soon as a 1 appears in Therefore, (7.131) holds. 

Estimating \gi\ < we see that, as in the proof of Lemma 7.13 that any nonzero term on the 
RHS of (7.129) satisfies 

■ ■ 9^2rQ^2r{F^,^)^2r\ ^ + (7.133) 


38 
























where s is the number of single labels appearing in the term. It then follows that, with / = |r|, 


2 l. ,lr 


fci=l fci=l 


< C 


-4r 

a -- 

_/Y2r ^ j^^r+s/2 




iV 2 r 


(7.134) 


□ 


where we used that s/2 + r — / > 0. As in the proof of Lemma 7.13 we now get the claim. 

7.3 Proof of Theorem 3.5 

In this section we show how the Helffer-Sjostrand formula is used to yield rigidity for the eigenvalues 
of H. For the most part the material in this section is standard; see, e.g., [8]. In this section we will 
always work on the event that the statement of Theorem 3.3 holds. Let / be a smooth compactly 
supported function. We have 


fW = I f 

^ JR 


iyf"(x)x(l/) + i(f(x) + wf'{x))x'{y) 

X — X — iy 


dxdy 


(7.135) 


where x is a smooth compactly supported function that is 1 in a neighborhood of the origin. Fix now 
E 2 e lEo,qG and let El := (recall the definition of By in Definition 2.1). Take / to satisfy 


where 


/(S) = 0, Ei[Ei-l,E2+vl f{E) = l, Ee[EyE2l 


77 = 




N 


(7.136) 


(7.137) 


We can take |/'(F1)|, |/"(F1)| < C for FI near Ei and |/'(F1)| < Crj and |/"(F1)| < Crj 2 for E near E 2 . 
For this smoothed out eigenvalue counting function we prove the following estimate. 

Lemma 7.19. Let f be as above. There exists a C 2 > 0 such that with {^,i')-high probability we have 


I f(X)dpE(X) - J /(A)dpfc(A) 


< 


F' 


C 2 $ 


N ’ 


(7.138) 


where denotes the empirical measure of Ht- 
Proof. Take x a smooth cut-off function to satisfy 

X(y) = l, \y\<N^^^y, x(l/) = 0, |y| > + 1. 


(7.139) 


It is no loss of generality to assume that By > 10. Let S = rriN — mic and p = pN — Pfc- From (7.135) 
we have 


J/(A)p(A)dA 


= Re 
< C 


;J/(A)p(A)dA 

Jy/"(a^)x(y)Im S{x + iy)dxdy 
+ C j\f{x)x'{y)\\lra S{x + iy)\ + \yf{x)x'{y)\\ReS{x + iy)|dxdy. (7.140) 
Since \S{x -I- iy)| < /{Ny) for y > 10 we immediately get 


I 


\y\\f{x)\\x\y)\^eS{x + iy)\dxdy < C 




,Cl? 


N 


(7.141) 
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where we used that \f'{x)\ < r]~^ is nonzero only on an interval of length r] near E 2 , and is bounded 
near Ei. Since x'iu) ^ 0 only for < \y\ < + 1 and f{x) 7 ^ 0 only on an interval of length 

2 n^^v have 

r n^Bv 1 

J l/(a;)x'(y)l|Im S{x + iy)|dxdy < (7.142) 

We bound the second line of (7.140) as 


2 //"(a;)x(y)Im S{x + iy)dxdy 

< 

y/"(a;)x(y)Im S{x + iy)dxdy 

J 


J\x-Ei\<2,\y\<10 

+ 

y/"(a;)x(y)Im -^(x + iy)dxdy 


J|a;-Ei|<2,|y|>10 

+ 

yf{x)x{y)^^ S{x + iy)dxdy 


J\x-E 2 \< 2 ,\y\<T] 

+ 

yf{x)x{y)^^ S{x + iy)dxdy 


J\x-E2\<2,\y\>ri 

= : 

A 1 +A 2 + A 3 + ^4 


We easily see that 

\Ai\<N-^^ (7.144) 

using that dist(Si, supp p) > . For A 2 we integrate by parts in x and then use d^Im S = —dyKeS 

to integrate by parts in y to obtain 


A 2 < 


f'ix)dy{yxiy))ReS{x + iy)dxdy 


\x-Ei\<2,\y\>10 

/'(x)10x(10)ReS'(x + 10i)dx 


'\x—Ei\<2 

Using I S'! < /{Ny) we see the first term is bounded by 

,log(A^)</?'^i^ 


= : Bi + B2. 


Bi < C- 


N 


(7.145) 


(7.146) 


The term B 2 is easily bounded by . 

For the term A^, note that since that y —> yim \m{E + \y)\ is a decreasing function for any Stieltjes 
transform m of a positive measure, we see that, using that is bounded near E 2 , 


Tf 

Im {rriMix + iy)] < 4.1m {rriMix + iy)] < C 

y 


Ny 


y <v 


(7.147) 


and so Im [S{x + iy)] < /{Ny) for y < Using then that the x integration is over an interval 

of length rj and \f"\ < (7 77 “^ we see that 


A 3 <Cr]<C 



For A 4 we integrate by parts as before and obtain 


VI 

f /'(x)dy(yx(y))Re5(x + iy)dxdy 


'^\y\>V,\^-B 2\<‘2 


+ 

/'(x)r 7 x( 7 )Re 5 (x + iy)dx 

— : B 3 + B 4 


J\x—E 2\<2 



(7.148) 


(7.149) 
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The term is estimated using l^l < {p^^^/{Ny). For the term we use that the x integration is 
over an interval of length rj, that \f'{x)\ < Crj~^ and |S| < ip'^^^/{Nr]) to obtain 


Bi<C 


~ir' 


(7.150) 


This yields the claim. □ 

Define the empirical eigenvalue counting function and the eigenvalue counting function of the 
deformed semicircle law by 

1 r® 

nt{E) = -\{i:Xi{Ht)<E}\, nfc,i(S) = J (7.151) 


respectively. From Lemma 7.19 we conclude the following lemma. 

Lemma 7.20. Let q > 0 be as above and let uj > 0 and t e Ti^. With [^,u)-high probability we have 
uniformly for E e lEo,qG, 

\nt{E) - nicAE)\ < (7.152) 

for a eonstant C 3 > 0 depending on Zy, Aq, Ci in (3.2), the constants appearing in (2.3) and the 
choice of q. 

Proof. Let / be as in Lemma 7.19. Note that spectrum of Ht is contained in — 1,N^^ + 1] 

with ((^, z/)-high probability. The left endpoint E^ of the support of pfc^t is the smallest solution to 
the equation 


T'l 


N {V^ -E_- Tmi,,T{E)A ■ 


(7.153) 


Clearly E- > — 1. With (^, z^)-high probability we then have that, with p as in Lemma 7.19, 

nt{E 2 ) < J f{E')dnt{E') < J f(E')dpf,AE') + ^ 

C2 ^ C2 ^ 

< nfc,t(E2 + p) + < nfc,t(E2) + + Cp 

< nfcAE 2 ) + ^ (7.154) 

where we used the fact that pfcg < C in ZEQ,qG- The lower bound is similar. □ 

The conclusion of Theorem 3.5 from Lemma 7.20 is standard and is similar to, for example, the 
proof of Theorem 7.6 in [ 8 ]. 


References 

[ 1 ] F. Bekerman, A. Figalli, and A. Guionnet. Transport maps for Beta-matrix models and univer¬ 
sality. preprint, arXiv:1311.2315, 2013. 

[2] P. Biane. On the free convolution with a semi-circular distribution. Indiana Univ. Math. J., 
46(3):705-718, 1997. 

[3] P. Bourgade, L. Erdos, and H.-T. Yau. Bulk universality of general /3-ensembles with non-convex 
potential. J. Math. Phys., 53(9);095221, 2012. 

[4] P. Bourgade, L. Erdos, and H.-T. Yau. Edge universality of /3-ensembles. Commun. Math. Phys., 
332(l);261-353, 2014. 

[5] P. Bourgade, L. Erdos, and H.-T. Yau. Universality of general /3-ensembles. Duke Math. J.l, 
163(6):1127--1190, 2014. 


41 









[6] P. Bourgade, L. Erdos, H.-T. Yau, and J. Yin. Fixed energy universality for generalized Wigner 
matrices, preprint, arXiv: 1407.5606, 2014. 

[7] L. Erdos, A. Knowles, H.-T. Yau, and J. Yin. Spectral statistics of Erdos-Renyi graphs II: 
Eigenvalue spacing and the extreme eigenvalues. Commun. in Math. Phys., 314(3):587-640, 
2012. 

[8] L. Erdos, A. Knowles, H.-T. Yau, and J. Yin. The local semicircle law for a general class of 
random matrices. Electron. J. Probah., 18(59): 1-58, 2013. 

[9] L. Erdos, A. Knowles, H.-T. Yau, and J. Yin. Spectral statistics of Erdos-Renyi graphs I: local 
semicircle law. Ann. Probah., 41(3B):2279-2375, 2013. 

[10] L. Erdos, S. Peche, J. A. Ramirez, B. Schlein, and H.-T. Yau. Bulk universality for wigner 
matrices. Comm. Pure Appl. Math., 63(7):895-925, 2010. 

[11] L. Erdos, B. Schlein, and H.-T. Yau. Wegner estimate and level repulsion for Wigner random 
matrices. Int. Math. Res. Not., 2010(3):436-479, 2010. 

[12] L. Erdos, B. Schlein, and H.-T. Yau. Universality of random matrices and local relaxation flow. 
Invent. Math., 185(1):75-119, 2011. 

[13] L. Erdos, B. Schlein, H.-T. Yau, and J. Yin. The local relaxation flow approach to universality 
of the local statistics for random matrices. Ann. Inst. Henri Poincare Probab. Stat., 48(l):l-46, 
2012. 

[14] L. Erdos and K. Schnelli. Universality for random matrix flows with time-dependent density. 
preprint, arxiv: 1504-00650. 

[15] L. Erdos and H.-T. Yau. Gap universality of generalized Wigner and /3-ensembles, preprint, 
arXiv: 1211.3786, 2012. 

[16] L. Erdos and H.-T. Yau. Universality of local spectral statistics of random matrices. Bull. Amer. 
Math. Soc., 49(3):377-414, 2012. 

[17] L. Erdos, H.-T. Yau, and J. Yin. Universality for generalized Wigner matrices with Bernoulli 
distribution. J. of Combinatorics, (2):15-85, 2011. 

[18] L. Erdos, H.-T. Yau, and J. Yin. Bulk universality for generalized Wigner matrices. Probab. 
Theory Related Fields, 154(l-2):341-407, 2012. 

[19] L. Erdos, H.-T. Yau, and J. Yin. Rigidity of eigenvalues of generalized Wigner matrices. Adv. 
Math., 229(3):1435-1515, 2012. 

[20] J. Huang, B. Landon, and H.-T. Yau. Bulk universality of sparse random matrices. J. Math. 
Phys., to appear, 2015. 

[21] J. O. Lee and K. Schnelli. Local deformed semicircle law and complete delocalization for Wigner 
matrices with random potential. J. Math. Phys., 54(10):103504, 2013. 

[22] J. O. Lee, K. Schnelli, B. Stetler, and H.-T. Yau. Bulk universality for deformed Wigner matrices. 
preprint, arXiv:1405.6634, 2014. 

[23] M. L. Mehta. Random matriees, volume 142. Academic press, 2004. 

[24] H. Nguyen, T. Tao, and V. Vu. Random matrices: tail bounds for gaps between eigenvalues. 
preprint, arXiv:1504:00396, 2015. 

[25] M. Shcherbina. Change of variables as a method to study general /3-models: bulk universality. J. 
Math. Phys., 55(4):043504, 2014. 


42 


[26] T. Tao and V. Vu. Random matrices: Universality of local eigenvalue statistics up to the edge. 
Commun. in Math. Phys., 298(2):549-572, 2010. 

[27] T. Tao and V. Vu. Random matrices: universality of local eigenvalue statistics. Acta Math., 
206(l):127-204, 2011. 

[28] E. Wigner. Characteristic vectors of bordered matrices with infinite dimensions. Ann. Math, 
62:548-564, 1955. 

[29] H.-T. Yau. Relative entropy and hydrodynamics of Ginzburg-Landau models. Lett. Math. Phys., 
22(l):63-80, 1991. 


43 


